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Physics invites the idea that space contains energy whose gravitational effect approximates that of 
Einstein's cosmological constant, A; nowadays the concept is termed dark energy or quintessence. 
Physics also suggests the dark energy could be dynamical, allowing the arguably appealing picture 
that the dark energy density is evolving to its natural value, zero, and is small now because the 
expanding universe is old. This alleviates the classical problem of the curious energy scale of order 
a millielectronvolt associated with a constant A. Dark energy may have been detected by recent 
advances in the cosmological tests. The tests establish a good scientific case for the context, in 
the relativistic Friedmann-Lemaitre model, including the gravitational inverse square law applied 
to the scales of cosmology. We have well-checked evidence that the mean mass density is not 
much more than one quarter of the critical Einstein-de Sitter value. The case for detection of dark 
energy is serious but not yet as convincing; we await more checks that may come out of work in 
progress. Planned observations might be capable of detecting evolution of the dark energy density; 
a positive result would be a considerable stimulus to attempts to understand the microphysics 
of dark energy. This review presents the basic physics and astronomy of the subject, reviews 
the history of ideas, assesses the state of the observational evidence, and comments on recent 
developments in the search for a fundamental theory. 
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I. INTRODUCTION 

There is significant observational evidence for the detection of Einstein's cosmological constant, A, or a component 
of the material content of the universe that varies only slowly with time and space and so acts like A. We will use 
the term dark energy for A or a component that acts like it. Detection of dark energy would be a new clue to an old 
puzzle, the gravitational effect of the zero-point energies of particles and fields. The total with other energies that are 
close to homogeneous and nearly independent of time act as dark energy. The puzzle has been that the value of the 
dark energy density has to be tiny compared to what is suggested by dimensional analysis; the startling new evidence 
is that it may be different from the only other natural value, zero. 

The main question to consider now has to be whether to accept the evidence for detection of dark energy. We 
outline the nature of the case in this section. After reviewing the basic concepts of the relativistic world model, in 
Sec. II, and in Sec. Ill reviewing the history of ideas, we present in Sec. IV a more detailed assessment of the state 
of the cosmological tests and the evidence for detection of A or its analog in dark energy. 

There is little new to report on the big issue for physics — why the dark energy density is so small — since 
Weinberg's (1989) review in this Journal 1 . But there have been analyses of a simpler idea: can we imagine the 
present dark energy density is evolving, perhaps approaching zero? Models are introduced in Sees. II. C and III.E, 
and recent work is summarized in more detail in the Appendix. Feasible advances in the cosmological tests could 
detect evolution of the dark energy density, and maybe its gravitational response to the large-scale fluctuations in the 
mass distribution. That would really drive the search for a more fundamental physics model for dark energy. 

A. The issues for observational cosmology 

We have to make two points. First, cosmology has a substantial observational and experimental basis that shows 
many aspects of the standard model almost certainly are good approximations to reality. Second, the empirical basis 
is not nearly as strong as it is for the standard model for particle physics: in cosmology it is not yet a matter of 
measuring the parameters in a well-established theory. 

To explain the second point we must remind those more accustomed to experiments in the laboratory than obser- 
vations in astronomy of the astronomers' tantalus principle: one can look at distant objects but never touch them. 
For example, the observations of supernovae in distant galaxies offer evidence for the detection of dark energy, under 
the assumption that distant and nearby supernovae are drawn from the same statistical sample (that is, that they 
are statistically similar enough for the purpose of this test). There is no direct way to check this, and it is easy 
to imagine differences between distant and nearby supernovae of the same nominal type. More distant supernovae 
are seen in younger galaxies, because of the light travel time, and these younger galaxies tend to have more massive 
rapidly evolving stars with lower heavy element abundances. How do we know the properties of the supernovae are 
not also different? We recommend Leibundgut's (2001, Sec. 4) discussion of the astrophysical hazards. Astronomers 
have checks for this and other issues of interpretation of the observations used in the cosmological tests. But it takes 
nothing away from this careful and elegant work to note that the checks seldom can be convincing, because the astron- 
omy is complicated and what can be observed is sparse. What is more, we don't know ahead of time that the physics 



1 Sahni and Starobinsky (2000), Carroll (2001), Witten (2001), Weinberg (2001), and Ellwanger (2002) present more recent reviews. 
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that is well tested on scales ranging from the laboratory to the Solar System survives the enormous extrapolation to 
cosmology. 

Our first point is that the situation is by no means hopeless, because we now have significant cross-checks from 
the consistency of results based on independent applications of the astronomy and of the physics of the cosmological 
model. If the physics or astronomy were faulty we would not expect consistency from independent lines of evidence 
— apart from the occasional accident, and the occasional tendency to stop the analysis when it approachs the "right 
answer". We have to demand abundant evidence of consistency, and that is starting to appear. 

The case for detection of A or dark energy commences with the Friedmann-Lemaitre cosmological model. In 
this model the expansion history of the universe is determined by a set of dimensionless parameters whose sum is 
normalized to unity, 

f^MO + ^RO + ^AO + ^KO = 1- (1) 

The first, SImo, is a measure of the present mean mass density in nonrelativistic matter, mainly baryons and nonbary- 
onic dark matter. The second, fiRo ~ 1 x 1CP 4 , is a measure of the present mass in the relativistic 3 K thermal cosmic 
microwave background radiation that almost homogeneously fills space, and the accompanying low mass neutrinos. 
The third is a measure of A or the present value of the dark energy equivalent. The fourth, £Iko, is an effect of the 
curvature of space. We review some details of these parameters in the next section, and of their measurements in 
Sec. IV. 

The most direct evidence for detection of dark energy comes from observations of supernovae of a type whose 
intrinsic luminosities are close to uniform (after subtle astronomical corrections, a few details of which are discussed 
in Sec. IV. B. 4). The observed brightness as a function of the wavelength shift of the radiation probes the geometry 
of spacetime, in what has come to be called the redshift-magnitude relation. 2 The measurements agree with the 
relativistic cosmological model with fixo = 0, meaning no space curvature, and Q\o ~ 0.7, meaning nonzero A. A 
model with fi\o = is two or three standard deviations off the best fit, depending on the data set and analysis 
technique. This is an important indication, but 2 to 3 a is not convincing, even when we can be sure the systematic 
errors are under reasonable control. And we have to consider that there may be a significant systematic error from 
differences between distant, high redshift, and nearby, low redshift, supernovae. 

There is a check, based on the CDM model 3 for structure formation. The fit of the model to the observations 
reviewed in Sec. IV. B yields two key constraints. First, the angular power spectrum of fluctuations in the temperature 
of the 3 K thermal cosmic microwave background radiation across the sky indicates Oko is small. Second, the power 
spectrum of the spatial distribution of the galaxies requires Qmo ~ 0.25. Similar estimates of Qmo follow from 
independent lines of observational evidence. The rate of gravitational lcnsing prefers a somewhat larger value (if f2xo 
is small), and some dynamical analyses of systems of galaxies prefer lower J7mo- But the differences could all be in 
the measurement uncertainties. Since J7ro in Eq. (|l]) is small, the conclusion is £Iao is large, in excellent agreement 
with what the supernovae say. 

Caution is in order, however, because this check depends on the CDM model for structure formation. We can't see 
the dark matter, so we naturally assign it the simplest properties we can get away with. Maybe it is significant that 
the model has observational problems with galaxy formation, as discussed in Sec. IV. A. 2, or maybe these problems 
are only apparent, from the complications of the astronomy. We are going to have to determine which it is before we 
can have a lot of confidence in the role of the CDM model in the cosmological tests. We will get a strong hint from 
the precision measurements in progress of the angular distribution of the 3 K thermal cosmic microwave background 
radiation. 4 If the results match in all detail the prediction of the relativistic model for cosmology and the CDM model 
for structure formation, with parameter choices that agree with the constraints from all the other cosmological tests, 
it will be strong evidence that we are approaching a good approximation to reality, and the completion of the great 
program of cosmological tests that commenced in the 1930s. But all that is to come. 

We emphasize that the advances in the empirical basis for cosmology already are very real and substantial. How 
firm is the conclusion depends on the issue, of course. Every competent cosmologist we know accepts as established 
beyond reasonable doubt that the universe is expanding and cooling in a near homogeneous and isotropic way from a 
hotter denser state: how else could space, which is transparent now, have been filled with radiation that has relaxed 
to a thermal spectrum? The debate is about when the expansion commenced or became a meaningful concept. Some 



2 The apparent magnitude is m = — 2.51og 10 / plus a constant, where / is the detected energy flux density in a chosen wavelength band. 
The standard measure of the wavelength shift, due to the expansion of the universe, is the redshift z defined in Eq. {^) below. 

3 The model is named after the nonbaryonic cold dark matter (CDM) that is assumed to dominate the masses of galaxies in the present 
universe. There are more assumptions in the CDM model, of course; they are discu ssed in Sees. TTT.D and TV. A .2. 

4 At the time of writing the MAP satellite is taking data; the project is described in fittp:/ /map. gsfc.nasa.gov/. 
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whose opinions and research we respect question the extrapolation of the gravitational inverse square law, in its use in 
estimates of masses in galaxies and systems of galaxies, and of Quo- We agree that this law is one of the hypotheses 
to be tested. Our conclusion from the cosmological tests in Sec. IV is that the law passes significant though not yet 
complete tests, and that we already have a strong scientific case, resting on the abundance of cross-checks, that the 
matter density parameter ilpjo is about one quarter. The case for detection of J} A o is significant too, but not yet as 
compelling. 

For the most part the results of the cosmological tests agree wonderfully well with accepted theory. But the obser- 
vational challenges to the tests are substantial: we are drawing profound conclusions from very limited information. 
We have to be liberal in considering ideas about what the universe is like, conservative in accepting ideas into the 
established canon. 



B. The opportunity for physics 

Unless there is some serious and quite unexpected flaw in our understanding of the principles of physics we can 
be sure the zero-point energy of the electromagnetic field at laboratory wavelengths is real and measurable, as in the 
Casimir (1948) effect. 5 Like all energy, this zero-point energy has to contribute to the source term in Einstein's grav- 
itational field equation. If, as seems likely, the zero-point energy of the electromagnetic field is close to homogeneous 
and independent of the velocity of the observer, it manifests itself as a positive contribution to Einstein's A, or dark 
energy. The zero point energies of the fermions make a negative contribution. Other contributions, perhaps including 
the energy densities of fields that interact only with themselves and gravity, might have either sign. The value of the 
sum suggested by dimensional analysis is much larger than what is allowed by the relativistic cosmological model. 
The only other natural value is A = 0. If A really is tiny but not zero, it adds a most stimulating though enigmatic 
clue to physics to be discovered. 

To illustrate the problem we outline an example of a contribution to A. The energy density in the 3 K thermal 
cosmic microwave background radiation, which amounts to f^Ro ~ 5 x 1CP 5 in Eq. (Q) (ignoring the neutrinos) peaks 
at wavelength A ~ 2 mm. At this Wien peak the photon occupation number is of order a fifteenth. The zero-point 
energy amounts to half the energy of a photon at the given frequency. This means the zero-point energy in the 
electromagnetic field at wavelengths A ~ 2 mm amounts to a contribution SQao ~4x 10~ 4 to the density parameter 
in A or dark energy. The sum over modes scales as A -4 (as illustrated in Eq. p7|). Thus a naive extrapolation to 
visible wavelengths says the contribution amounts to Sil\o ~ 5 x 10 10 , already a ridiculous number. 

The situation can be compared to the development of the theory of the weak interactions. The Fermi point-like 
interaction model is strikingly successful for a considerable range of energies, but it was clear from the start that 
the model fails at high energy. A fix was discussed — mediate the interaction by an intermediate boson — and 
eventually incorporated into the even more successful electroweak theory. General relativity and quantum mechanics 
are strikingly successful on a considerable range of length scales, provided we agree not to use the rules of quantum 
mechanics to count the zero-point energy density in the vacuum, even though we know we have to count the zero-point 
energies in all other situations. There are thoughts on how to improve the situation, though they seem less focused 
than was the case for the Fermi model. Maybe a new energy component spontaneously cancels the vacuum energy 
density; maybe the new component varies slowly with position and here and there happens to cancel the vacuum 
energy density well enough to allow observers like us to flourish. Whatever the nature of the more perfect theory, it 
must reproduce the successes of general relativity and quantum mechanics. That includes the method of representing 
the material content of the observable universe — all forms of mass and energy — by the stress-energy tensor, and the 
relation between the stress-energy tensor and the curvature of macroscopic spacetime. One part has to be adjusted. 

The numerical values of the parameters in Eq. ([!]) also are enigmatic, and possibly trying to tell us something. The 
evidence is that the parameters have the approximate values 

O\ ~0.7, rbMo~0.2, ft BO ~0.05. (2) 



5 See Bordag, Mohideen, and Mostepanenko (2001) for a recent review. The attractive Casimir force between two parallel conducting 
plates results from the boundary condition that suppresses the number of modes of oscillation of the electromagnetic field between the 
plates, thus suppressing the energy of the system. One can understand the effect at small separation without reference to the quantum 
behavior of the electromagnetic field, as in the analysis of the Van der Waals interaction in quantum mechanics, by taking account of the 
term in the particle Hamiltonian for the Coulomb potential energy between the charged particles in the separate neutral objects. But a 
more complete treatment, as discussed by Cohen- Tannoudji, Dupont-Roc, and Grynberg (1992), replaces the Coulomb interaction with 
the coupling of the charged particles to the electromagnetic field operator. In this picture the Van der Waals interaction is mediated by 
the exchange of virtual photons. In either way of looking at the Casimir effect — the perturbation of the normal modes or the exchange 
of virtual quanta of the unperturbed modes — the effect is the same, the suppression of the energy of the system. 
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We have written JImo i n two parts: Qbo measures the density of the baryons we know exist and Odmo that of the 
hypothetical nonbaryonic cold dark matter we need to fit the cosmological tests. The parameters fieo an d ^dmo have 
similar values but represent different things — baryonic and nonbaryonic matter — and f^Ao, which is thought to 
represent something completely different, is not much larger. Also, if the parameters were measured when the universe 
was one tenth its present size the time-independent A parameter would contribute J7a ~ 0.003. That is, we seem to 
have come on the scene just as A has become an important factor in the expansion rate. These curiosities surely are 
in part accidental, but maybe in part physically significant. In particular, one might imagine that the dark energy 
density represented by A is rolling to its natural value, zero, and is very small now because we measure it when the 
universe is very old. We will discuss efforts along this line to at least partially rationalize the situation. 

C. Some explanations 

We have to explain our choice of nomenclature. Basic concepts of physics say space contains homogeneous zero- 
point energy, and maybe also energy that is homogeneous or nearly so in other forms, real or effective (as from counter 
terms in the gravity physics, which make the net energy density cosmologically acceptable). In the literature this 
near homogeneous energy has been termed the vacuum energy, the sum of vacuum energy and quintessence (Caldwell, 
Dave, and Steinhardt, 1998), and the dark energy (Turner, 1999). We have adopted the last term, and we will refer 
to the dark energy density pa that manifests itself as an effective version of Einstein's cosmological constant, but one 
that may vary slowly with time and position. 6 

Our subject involves two quite different traditions, in physics and astronomy. Each has familiar notation, and 
familiar ideas that may be "in the air" but not in the recent literature. Our attempt to take account of these 
traditions commences with the summary in Sec. II of the basic notation with brief explanations. We expect readers 
will find some of these concepts trivial and others of some use, and that the useful parts will be different for different 
readers. 

We offer in Sec. Ill our reading of the history of ideas on A and its generalization to dark energy. This is a 
fascinating and we think edifying illustration of how science may advance in unexpected directions. It is relevant to 
an understanding of the present state of research in cosmology, because traditions inform opinions, and people have 
had mixed feelings about A ever since Einstein (1917) introduced it 85 years ago. The concept never entirely dropped 
out of sight in cosmology because a series of observations hinted at its presence, and because to some cosmologists 
A fits the formalism too well to be ignored. The search for the physics of the vacuum, and its possible relation to 
A, has a long history too. Despite the common and strong suspicion that A must be negligibly small, because any 
other acceptable value is absurd, all this history has made the community well prepared for the recent observational 
developments that argue for the detection of A. 

Our approach in Sec. IV to the discussion of the evidence for detection of A, from the cosmological tests, also 
requires explanation. One occasionally reads that the tests will show us how the world ends. That certainly seems 
interesting, but it is not the main point: why should we trust an extrapolation into the indefinite future of a theory 
we can at best show is a good approximation to reality? 7 As we remarked in Sec. LA, the purpose of the tests is to 
check the approximation to reality, by checking the physics and astronomy of the standard relativistic cosmological 
model, along with any viable alternatives that may be discovered. We take our task to be to identify the aspects of 
the standard theory that enter the interpretation of the measurements and thus are or may be empirically checked or 
measured. 

II. BASIC CONCEPTS 

A. The Friedmann-LemaTtre model 

The standard world model is close to homogeneous and isotropic on large scales, and lumpy on small scales — the 
effect of the mass concentrations in galaxies, stars, people, and all that. The length scale at the transition from nearly 
smooth to strongly clumpy is about 10 Mpc. We use here and throughout the standard astronomers' length unit, 

1 Mpc = 3.1 x 10 24 cm = 3.3 x 10 6 light years. (3) 



6 The dark energy should of course be distinguished from a hypothetical gas of particles with velocity dispersion large enough that the 
distribution is close to homogeneous. 

7 Observations may now have detected A, at a characteristic energy scale of a millielectronvolt (Eq. |47| ) . We have no guarantee that 
there is not an even lower energy scale; such a scale could first become apparent through the cosmological tests. 
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To be more definite, imagine many spheres of radius 10 Mpc are placed at random, and the mass within each is 
measured. At this radius the rms fluctuation in the set of values of masses is about equal to the mean value. 
On smaller scales the departures from homogeneity are progressively more nonlinear; on larger scales the density 
fluctuations are perturbations to the homogeneous model. From now on we mention these perturbations only when 
relevant for the cosmological tests. 

The expansion of the universe means the distance l(t) between two well-separated galaxies varies with world time, 
t, as 

l(t) (x a(t), (4) 

where the expansion or scale factor, a(t), is independent of the choice of galaxies. It is an interesting exercise, for 
those who have not already thought about it, to check that Eq. (||) is required to preserve homogeneity and isotropy. 8 
The rate of change of the distance in Eq. (||) is the speed 

v = dl/dt = HI, H = a/a, (5) 

where the dot means the derivative with respect to world time t and H is the time-dependent Hubble parameter. 
When v is small compared to the speed of light this is Hubble's law. The present value of H is Hubble's constant, 
Hq. When needed we will use 9 

Ho = 100ft km s^Mpc -1 = 67 ± 7 km s^Mpc -1 = (15 ± 2 Gyr)" 1 , (6) 

at two standard deviations. The first equation defines the dimensionless parameter h. 

Another measure of the expansion follows by considering the stretching of the wavelength of light received from 
a distant galaxy. The observed wavelength, A b s , of a feature in the spectrum that had wavelength A em at emission 
satisfies 



l + z = A obs /A om = a(t ohs )/a(t cm ), (7) 

where the expansion factor a is defined in Eq. (Q) and z is the redshift. That is, the wavelength of freely traveling 
radiation stretches in proportion to the factor by which the universe expands. To understand this, imagine a large 
part of the universe is enclosed in a cavity with perfectly reflecting walls. The cavity expands with the general 
expansion, the widths proportional to a(t). Electromagnetic radiation is a sum of the normal modes that fit the 
cavity. At interesting wavelengths the mode frequencies are much larger than the rate of expansion of the universe, so 
adiabaticity says a photon in a mode stays there, and its wavelength thus must vary as A cx a(t), as stated in Eq. (Q). 
The cavity disturbs the long wavelength part of the radiation, but the disturbance can be made exceedingly small by 
choosing a large cavity. 

Equation defines the redshift z. The redshift is a convenient label for epochs in the early universe, where z 
exceeds unity. A good exercise for the student is to check that when z is small Eq. (|7]) reduces to Hubble's law, where 
Xz is the first-order Doppler shift in the wavelength A, and Hubble's parameter H is given by Eq. (^). Thus Hubble's 
law may be written as cz = HI (where we have put in the speed of light). 

These results follow from the symmetry of the cosmological model and conventional local physics; wc do not need 
general relativity theory. When z ^ 1 we need the relativistic theory to compute the relations among the redshift and 
other observables. An example is the relation between redshift and apparent magnitude used in the supernova test. 
Other cosmological tests check consistency among these relations, and this checks the world model. 



8 We feel we have to comment on a few details about Eq. (^) to avoid contributing to debates that are more intense than seem warranted. 
Think of the world time t as the proper time kept by each of a dense set of observers, each moving so all the others are isotropically 
moving away, and with the times synchronized to a common energy density, pit), in the near homogeneous expanding universe. The 
distance l(t) is the sum of the proper distances between neighboring observers, all measured at time t, and along the shortest distance 
between the two observers. The rate of increase of the distance, dl/dt, may exceed the velocity of light. This is no more problematic in 
relativity theory than is the large speed at which the beam of a flashlight on Earth may swing across the face of the Moon (assuming an 
adequately tight beam). Space sections at fixed t may be non-compact, and the total mass of a homogeneous universe formally infinite. 
As far as is known this is not meaningful: we can only assert that the universe is close to homogeneous and isotropic over observable 
scales, and that what can be observed is a finite number of baryons and photons. 

9 The numerical values in Eq. fet) are determined from an analysis of almost all published measurements of Ho prior to mid 1999 (Gott 
et al., 2001). They are a very reasonable summary of the current situation. For instance, the Hubble Space Telescope Key Project 
summary measurement value Hq = 72 ± 8 km s _1 Mpc _1 (1 a uncertainty, Freedman et al., 2001) is in very good agreement with 
Eq. (Ej), as is the recent Tammann et al. (2001) summary value Ho = 60 ± 6 km s — 1 Mpc~ 1 (approximate 1 a systematic uncertainty). 
This is an example of the striking change in the observational situation over the previous 5 years: the uncertainty in Ho has decreased 
by more than a factor of 3, making it one of the better measured cosmological parameters. 
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In general relativity the second time derivative of the expansion factor satisfies 

- = -^G(p + 3p). (8) 
a 3 

The gravitational constant is G. Here and throughout we choose units to set the velocity of light to unity. The mean 
mass density, p(t), and the pressure, p(t), counting all contributions including dark energy, satisfy the local energy 
conservation law, 

p = -3~(p+p). (9) 

The first term on the right-hand side represents the decrease of mass density due to the expansion that more broadly 
disperses the matter. The pdV work in the second term is a familiar local concept, and meaningful in general relativity. 
But one should note that energy does not have a general global meaning in this theory. 
The first integral of Eqs. (||) and (^) is the Friedmann equation 

g 

a 2 = -irGpa 2 + constant. (10) 
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It is conventional to rewrite this as 

2 



Q =HiE(z) 2 = H 2 (n MO (i + z) 3 + Q RO (i + z) 4 + n A o + n KO {i + z) 2 ). (11) 

The first equation defines the function E(z) that is introduced for later use. The second equation assumes constant 
A; the time-dependent dark energy case is reviewed in Sees. II. C and III.E. The first term in the last part of Eq. ( |Tl| ) 
represents non-relativistic matter with negligibly small pressure; one sees from Eqs. (@) and (|) that the mass density 
in this form varies with the expansion of the universe as pm oc a~ 3 oc (1 + z) 3 . The second term represents radiation 
and relativistic matter, with pressure = pr/3, whence oc (1 + z) 4 . The third term is the effect of Einstein's 
cosmological constant, or a constant dark energy density. The last term, discussed in more detail below, is the constant 
of integration in Eq. (|l0|). The four density parameters fiio are the fractional contributions to the square of Hubble's 
constant, Hq, that is, flio(t) — 8nG pio / '(3Hq) . At the present epoch, z — 0, the present value of a/a is Hq, and the 
fliO sum to unity (Eq. p|) . 
In this notation, Eq. (H) is 



a 
a 



-h 2 (n M0 (i + zf/2 + n K0 (i + zf - n A0 ) . (12) 



The constant of integration in Eqs. (10) and (jllj) is related to the geometry of spatial sections at constant world 
time. Recall that in general relativity events in spacetime are labeled by the four coordinates x^ of time and space. 
Neighboring events 1 and 2 at separation dx* 1 have invariant separation ds defined by the line element 

ds 2 = g^dx^dx". (13) 

The repeated indices are summed, and the metric tensor is a function of position in spacetime. If ds 2 is positive 
then ds is the proper (physical) time measured by an observer who moves from event 1 to 2; if negative, \ds\ is 
the proper distance between events 1 and 2 measured by an observer who is moving so the events are seen to be 
simultaneous. 

In the flat spacetime of special relativity one can choose coordinates so the metric tensor has the Minkowskian form 












-1 











-1 











-1 



V. = o o _! o • (14) 
\0 -1/ 

A freely falling, inertial, observer can choose locally Minkowskian coordinates, such that along the path of the observer 
9nv = f]{iv and the first derivatives of c/ M „ vanish. 

In the homogeneous world model we can choose coordinates so the metric tensor is of the form that results in the 
line element 

ds 2 - dt 2 - a(t) 2 ( ^ „ + r 2 (dB 2 + sin 2 

V i + aH 

= dt 2 - K- X a{tf (d X 2 + sinh 2 x(d0 2 + sin 2 Odcf 2 )) . (15) 
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In the second expression, which assumes K > 0, the radial coordinate is r = K^ 1 / 2 sinhx. The expansion factor 
a(t) appears in Eq. (Q). If a were constant and the constant K vanished this would represent the flat spacetime of 
special relativity in polar coordinates. The key point for now is that SIko m Eq. (]ll]), which represents the constant 
of integration in Eq. ([To|), is related to the constant K: 

n KQ = K/(H a a ) 2 , (16) 

where ao is the present value of the expansion factor a(t). Cosmological tests that are sensitive to the geometry of 
space constrain the value of the parameter J^koi and ^ko and the other density parameters flio in Eq. ( pd| ) determine 
the expansion history of the universe. 

It is useful for what follows to recall that the metric tensor in Eq. satisfies Einstein's field equation, a differential 
equation we can write as 

G>„ — SnGT^. (17) 

The left side is a function of g^ v and its first two derivatives and represents the geometry of spacetime. The stress- 
energy tensor T^ v represents the material contents of the universe, including particles, radiation, fields, and zero-point 
energies. An observer in a homogeneous and isotropic universe, moving so the universe is observed to be isotropic, 
would measure the stress-energy tensor to be 
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(18) 



This diagonal form is a consequence of the symmetry; the diagonal components define the pressure and energy density. 
With Eq. (|l8|), the differential equation ( |l7j ) yields the expansion rate equations ( pd| ) and (p^). 



B. The cosmological constant 

Special relativity is very successful in laboratory physics. Thus one might guess any inertial observer would see the 
same vacuum. A freely moving inertial observer represents spacetime in the neighborhood by locally Minkowskian 
coordinates, with the metric tensor rj^y given in Eq. (|l4|). A Lorentz transformation to an inertial observer with 
another velocity does not change this Minkowski form. The same must be true of the stress-energy tensor of the 
vacuum, if all observers see the same vacuum, so it has to be of the form 

T^v = PA9^, (19) 

where p\ is a constant, in a general coordinate labeling. On writing this contribution to the stress-energy tensor 
separately from all the rest, we bring the field equation (jlj]) to 

G^^8ttG(T^+ P a9^). (20) 

This is Einstein's (1917) revision of the field equation of general relativity, where pa is proportional to his cosmological 
constant A; his reason for writing down this equation is discussed in Sec. III. A. In many dark energy scenarios p\ is 
a slowly varying function of time and its stress-energy tensor differs slightly from Eq. (|l^) , so the observed properties 
of the vacuum do depend on the observer's velocity. 

One sees from Eqs. (jl^), (|l8|), and (|l9|) that the new component in the stress-energy tensor looks like an ideal fluid 
with negative pressure 

Pa = -pA- (21) 
This fluid picture is of limited use, but the following properties are worth noting. 10 



These arguments have been familiar, in some circles, for a long time, though in our experience discussed more often in private than the 
literature. Early statements of elements are in Lemaitre (1934) and McCrea (1951); see Kragh (1999, pp. 397-8) for a brief historical 
account. 
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The stress-energy tensor of an ideal fluid with four-velocity u M generalizes from Eq. (|lj) to T^ u = (p+p)u^u v —pg^ . 
The equations of fluid dynamics follow from the vanishing of the divergence of T M1/ . Let us consider the simple case of 
locally Minkowskian coordinates, meaning free fall, and a fluid that is close to homogeneous. By the latter we mean 
the fluid velocity v — the space part of the four-velocity — and the density fluctuation dp from homogeneity may 
be treated in linear perturbation theory. Then the equations of energy and momentum conservation are 

5p+((p) + (p))V -v = 0, ((p) + ( P ))j+c 2 s V6p = Q, (22) 

where c 2 = dp /dp and the mean density and pressure are (p) and (p). These combine to 

Sp = c 2 V 2 Sp. (23) 

If c 2 is positive this is a wave equation, and c s is the speed of sound. 

The first of Eqs. ( p2| ) is the local energy conservation law, as in Eq. (||). If p = —p, the pdV work cancels the 
pdV part: the work done to increase the volume cancels the effect of the increased volume. This has to be so for a 
Lorentz-invariant stress-energy tensor, of course, where all inertial observers see the same vacuum. Another way to 
see this is to note that the energy flux density in Eqs. (|22] ) is ((p) + (p))v. This vanishes when p = —p: the streaming 
velocity loses meaning. When c 2 is negative Eq. ( p3| ) says the fluid is unstable, in general. But when p — —p the 
vanishing divergence of T^ v becomes the condition seen in Eq. (|2^) that p = (p) + Sp is constant. 

There are two measures of gravitational interactions with a fluid: the passive gravitational mass density determines 
how the fluid streaming velocity is affected by an applied gravitational field, and the active gravitational mass density 
determines the gravitational field produced by the fluid. When the fluid velocity is nonrelativistic the expression for 
the former in general relativity is p+p, as one sees by writing out the covariant divergence of T^ v . This vanishes when 
p = —p, consistent with the loss of meaning of the streaming velocity. The latter is p + 3p, as one sees from Eq. (j^). 
Thus a fluid with p = — p/'i, if somehow kept homogeneous and static, would produce no gravitational field. 11 In the 
model in Eqs. ( |l9|) and ( pl|) the active gravitational mass density is negative when pa is positive. When this positive 
Pa dominates the stress-energy tensor a is positive: the rate of expansion of the universe increases. In the language 
of Eq. (^p| ) , this cosmic repulsion is a gravitational effect of the negative active gravitational mass density, not a new 
force law. 

The homogeneous active mass represented by A changes the equation of relative motion of freely moving test 
particles in the nonrelativistic limit to 



~d¥ 



n A0 H£r, (24) 



where g is the relative gravitational acceleration produced by the distribution of ordinary matter. 12 For an illustration 
of the size of the last term consider its effect on our motion in a nearly circular orbit around the center of the Milky 
Way galaxy. The Solar System is moving at speed v c = 220 km s _1 at radius r = 8 kpc. The ratio of the acceleration 
<7a produced by A to the total gravitational acceleration g — v 2 /r is 

gA/g = nAoH 2 r 2 /v 2 ~W- 5 , (25) 

a small number. Since we are towards the edge of the luminous part of our galaxy, a search for the effect of A on the 
internal dynamics of galaxies like the Milky Way does not look promising. The precision of celestial dynamics in the 
Solar System is much greater, but the effect of A is very much smaller; for the orbit of the Earth, 5a/.9 ~ 10~ 22 . 



One can generalize Eq. (19) to a variable pa, by taking pa to be negative but different from — pa- But if the 
dynamics were that of a fluid, with pressure a function of pa, stability would require c 2 = dp a /dp a > 0, from 
Eq. (^3|), which seems quite contrived. A viable working model for a dynamical pa is the dark energy of a scalar field 
with self-interaction potential chosen to make the variation of the field energy acceptably slow, as discussed next. 



Lest we contribute to a wrong problem for the student we note that a fluid with p = — p/3 held in a container would have net positive 
gravitational mass, from the pressure in the container walls required for support against the negative pressure of the contents. We have 
finessed the walls by considering a homogeneous situation. We believe Whittaker (1935) gives the first derivation of the relativistic 
active gravitational mass density. Whittaker also presents an example of the general proposition that the active gravitational mass of 
an isolated stable object is the integral of the time-time part of the stress-energy tensor in the locally Minkowskian rest frame. Misner 
and Putman (1959) give the general demonstration. 

This assumes the particles are close enough for application of the ordinary operational definition of proper relative position. The 
parameters in the last term follow from Eqs. (ph and (Ell). 
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C. Inflation and dark energy 

The negative active gravitational mass density associated with a positive cosmological constant is an early precursor 
of the inflation picture of the early universe; inflation in turn is one precursor of the idea that A might generalize into 
evolving dark energy. 

To begin, we review some aspects of causal relations between events in spacetime. Neglecting space curvature, a 
light ray moves proper distance dl = a(t)dx = dt in time interval dt, so the integrated coordinate displacement is 

dt/a(t). (26) 

If f^Ao = this integral converges in the past — we see distant galaxies that at the time of observation cannot have 
seen us since the singular start of expansion at a = 0. This "particle horizon problem" is curious: how could distant 
galaxies in different directions in the sky know to look so similar? The inflation idea is that in the early universe the 
expansion history approximates that of de Sitter's (1917) solution to Einstein's field equation for A > and T^ v = in 
Eq. (^p|). We can choose the coordinate labels in this de Sitter spacetime so space curvature vanishes. Then Eqs. ( |ll| ) 
and ( J12j ) say the expansion parameter is 

a cx e HKt , (27) 



where Ha is a constant. As one sees by working the integral in Eq. fl26|), here everyone can have seen everyone else in 
the past. The details need not concern us; for the following discussion two concepts are important. First, the early 
universe acts like an approximation to de Sitter's solution because it is dominated by a large effective cosmological 
"constant", or dark energy density. Second, the dark energy is modeled as that of a near homogeneous field, 

In this scalar field model, motivated by grand unified models of very high energy particle physics, the action of the 
real scalar field, $ (in units chosen so Planck's constant h is unity) is 



d 4 x. 



(28) 



The potential energy density V is a function of the field and g is the determinant of the metric tensor. When 
the field is spatially homogeneous (in the line element of Eq. [[l5|), and space curvature may be neglected, the field 
equation is 

* + 3 "* + S = °- ( 29 ) 
a d$ 

The stress-energy tensor of this homogeneous field is diagonal (in the rest frame of an observer moving so the universe 
is seen to be isotropic), with time and space parts along the diagonal 

p* = ^ 2 + y($), 

P<s> = \& (30) 

If the scalar field varies slowly in time, so that <I> 2 <C V, the field energy approximates the effect of Einstein's 
cosmological constant, with p$ ~ — 

The inflation picture assumes the near exponential expansion of Eq. ( p7| ) in the early universe lasts long enough 
that every bit of the present observable universe has seen every other bit, and presumably has discovered how to 
relax to almost exact homogeneity. The field $ may then start varying rapidly enough to produce the entropy of our 
universe, and the field or the entropy may produce the baryons, leaving p$ small or zero. But one can imagine the 
late time evolution of p$ is slow. If slower than the evolution in the mass density in matter, there comes a time when 
/9$ again dominates, and the universe appears to have a cosmological constant. 

A model for this late time evolution assumes a potential of the form 

V = k/$ c \ (31) 

where the constant k has dimensions of mass raised to the power a + 4. For simplicity let us suppose the universe 
after inflation but at high redshift is dominated by matter or radiation, with mass density p, that drives power law 
expansion, a cx t n . Then the power law solution to the field equation ( |29| ) with the potential in Eq. ( |3l| ) is 

$oci 2 /( 2+a \ (32) 
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and the ratio of the mass densities in the scalar field and in matter or radiation is 

p^/poct^ 2+a \ (33) 

In the limit where the parameter a approaches zero, p$ is constant, and this model is equivalent to Einstein's A. 

When a > the field $ in this model grows arbitrarily large at large time, so p$ — > 0, and the universe approaches 
the Minkowskian spacetime of special relativity. This is within a simple model, of course. It is easy to imagine that in 
other models p$ approaches a constant positive value at large time, and spacetime approaches the de Sitter solution, 
or p$ passes through zero and becomes negative, causing spacetime to collapse to a Big Crunch. 

The power law model with a > has two properties that seem desirable. First, the solution in Eq. (|32| ) is said to be 
an attractor (Ratra and Peebles, 1988) or tracker (Steinhardt, Wang, and Zlatev, 1999), meaning it is the asymptotic 
solution for a broad range of initial conditions at high redshift. That includes relaxation to a near homogeneous 
energy distribution even when gravity has collected the other matter into nonrelativistic clumps. Second, the energy 
density in the attractor solution decreases less rapidly than that of matter and radiation. This allows us to realize the 
scenario: after inflation but at high redshift the field energy density p$ is small so it does not disturb the standard 
model for the origin of the light elements, but eventually p$ dominates and the universe acts as if it had a cosmological 
constant, but one that varies slowly with position and time. We comment on details of this model in Sec III.E. 



III. HISTORICAL REMARKS 

These comments on what people were thinking are gleaned from the literature and supplemented by private dis- 
cussions and our own recollections. More is required for a satisfactory history of the subject, of course, but we hope 
we have captured the main themes of the discussion and the way the themes have evolved to the present appreciation 
of the situation. 



A. Einstein's thoughts 

Einstein disliked the idea of an island universe in asymptotically flat spacetime, because a particle could leave the 
island and move arbitrarily far from all the other matter in the universe, yet preserve all its inertial properties, which 
he considered a violation of Mach's idea of the relativity of inertia. Einstein's (1917) cosmological model accordingly 
assumes the universe is homogeneous and isotropic, on average, thus removing the possibility of arbitrarily isolated 
particles. Einstein had no empirical support for this assumption, yet it agrees with modern precision tests. There is 
no agreement on whether this is more than a lucky guess. 

Motivated by the observed low velocities of the then known stars, Einstein assumed that the large-scale structure 
of the universe is static. He introduced the cosmological constant to reconcile this picture with his general relativity 
theory. In the notation of Eq. (|l2|), one sees that a positive value of f^Ao can balance the positive values of flyio 
and Oro for consistency with a = 0. The balance is unstable: a small perturbation to the mean mass density or the 
mass distribution causes expansion or contraction of the whole or parts of the universe. One sees this in Eq. j24|): 
the mass distribution can be chosen so the two terms on the right hand side cancel, but the balance can be upset by 
redistributing the mass. 13 

Einstein did not consider the cosmological constant to be part of the stress-energy term: his form for the field 
equation (in the streamlined notation of Eq. p7|] ) is 

- SirGpAg^ = SirGT^. (34) 

The left hand side contains the metric tensor and its derivatives; a new constant of nature, A, appears in the addition 
to Einstein's original field equation. One can equally well put Einstein's new term on the right hand side of the 
equation, as in Eq. (p0|), and count pi^g^v as part of the source term in the stress-energy tensor. The distinction 



To help motivate the introduction of A, Einstein (1917) mentioned a modification of Newtonian gravity physics that might make the 
theory well defined when the mass distribution is homogeneous. In Einstein's example, similar to what was considered by Seeliger 
and Neumann in the mid 1890s, the modified field equation for the gravitational potential tp is V 2 (/3 — Ai/j = AttGpm- This allows the 
nonsingular homogeneous static solution ip = —4itGpm/^- In this example the potential for an isolated point mass is the Yukawa form, 
tp oc e~^ r Jr. Trautman (1965) points out that this is not the nonrelativistic limit of general relativity with the cosmological term. 
Rather, Eq. (E4) follows from V 2 ^ = AnG(pM — 2/9a), where the active gravitational mass density of the A term is pj\ + 3pA = —%PA- 
Norton (1999)rcviews the history of ideas of this Seeliger- Neumann Yukawa-type potential in gravity physics. 
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becomes interesting when pA takes part in the dynamics, and the field equation is properly written with p\, or its 
generalization, as part of the stress-energy tensor. One would then say that the differential equation of gravity physics 
has not been changed from Einstein's original form; instead there is a new component in the content of the universe. 

Having taken it that the universe is static, Einstein did not write down the differential equation for a(t), and so 
did not see the instability. Friedmann (1922, 1924) found the evolving homogeneous solution, but had the misfortune 
to do it before the astronomy became suggestive. Slipher's measurements of the spectra of the spiral nebulae - 
galaxies of stars — showed most are shifted toward the red, and Eddington (1924, pp. 161-2) remarked that that 
might be a manifestation of the second, repulsive, term in Eq. (p4|). Lemaitre (1927) introduced the relation between 
Slipher's redshifts and a homogeneous matter-filled expanding relativistic world model. He may have been influenced 
by Hubble's work that was leading to the publication in 1929 of the linear redshift-distance relation (eq. ||): as a 
graduate student at MIT he attended a lecture by Hubble. 

In Lemaitre's (1927) solution the expanding universe traces asymptotically back to Einstein's static case. Lemaitre 
then turned to what he called the primeval atom, and we would term a Big Bang model. This solution expands 
from densities so large as to require some sort of quantum treatment, passes through a quasi-static approximation 
to Einstein's solution, and then continues expanding to de Sitter's (1917) empty space solution. To modern tastes 
this "loitering" model requires incredibly special initial conditions, as will be discussed. Lemaitre liked it because the 
loitering epoch allows the expansion time to be acceptably long for Hubble's (1929) estimate of Hq, which is an order 
of magnitude high. 

The record shows Einstein never liked the A term. His view of how general relativity might fit Mach's principle 
was disturbed by de Sitter's (1917) solution to Eq. @ for empty space (T^ = 0) with A > 0. 14 ' 15 Pais (1982, p. 
288) points out that Einstein in a letter to Weyl in 1923 comments on the effect of A in Eq. (|5|): "According to 
De Sitter two material points that are sufficiently far apart, continue to be accelerated and move apart. If there is 
no quasistatic world, then away with the cosmological term." We do not know whether at this time Einstein was 
influenced by Slipher's redshifts or Friedmann's expanding world model. 

The earliest published comments we have found on Einstein's opinion of A within the evolving world model (Einstein, 
1931; Einstein and de Sitter, 1932) make the point that, since not all the terms in the expansion rate Eq. (|j"l| ) are 
logically required, and the matter term surely is present and likely dominates over radiation at low redshift, a 
reasonable working model drops SIko and Q ao and ignores Q^q. This simplifies the expansion rate equation to what 
has come to be called the Einstein-de Sitter model, 

% = |ttGpm, (35) 

where pm is the mass density in non-relativistic matter; here 17m = 8itGpm/ (3H 2 ) is unity. The left side is a measure 
of the kinetic energy of expansion per unit mass, and the right-hand side a measure of the negative of the gravitational 
potential energy. In effect, this model universe is expanding with escape velocity. 

Einstein and de Sitter point out that Hubble's estimate of Hq and de Sitter's estimate of the mean mass density in 
galaxies are not inconsistent with Eq. (|3^) (and since both quantities scale with distance in the same way, this result 



North (1965) reviews the confused early history of ideas on the possible astronomical significance of de Sitter's solution for an empty 
universe with A > 0; we add a few comments on the physics that contributed to the discovery of the expanding world model. Suppose 
an observer in de Sitter's spacetime holds a string tied to a source of light, so the source stays at fixed physical distance r -C fTTJ . 
The source is much less massive than the observer, the gravitational frequency shift due to the observer's mass may be neglected, and 
the observer is moving freely. Then the observer receives light from the source shifted to the red by 5A/A = — (H\r) 2 /2. The observed 
redshifts of particles moving on geodesies depend on the initial conditions. Stars in the outskirts of our galaxy are held at fixed mean 
distances from us by their transverse motions. The mean shifts of the spectra of light from these stars include this quadratic de Sitter 
term as well as the much larger Doppler and ordinary gravitational shifts. The prescription for initial conditions that reproduces the 
linear redshift-distance relation for distant galaxies follows Weyl's (1923) principle: the world particle geodesies trace back to a near 
common position in the remote past, in the limiting case of the Friedmann-Lcmaitrc model at Qmo ~ ¥ 0- This spatially homogeneous 
coordinate labeling of de Sitter's spacetime, with space sections with negative curvature, already appears in de Sitter (1917, Eq. [15]), 
and is repeated in Lanczos (1922). This line element is the second expression in our Eq. ( ^j[ ) with a oc cosh Hj^t. Lemaitre (1925) and 
Robertson (1928) present the coordinate labeling for the spatially-flat case, where the line element is ds 2 = dt 2 — e 2HAt (dx 2 + dy 2 -\-dz 2 ) 
(in the choice of symbols and signature in Eqs. ||15j and PjJ). Lemaitre (1925) and Robertson (1928) note that particles at rest in this 
coordinate system present a linear redshift-distance relation, v = H^r, at small v. Robertson (1928) estimated H/±, and Lemaitre (1927) 
its analog for the Friedmann-Lemaitre model, from published redshifts and Hubble's galaxy distances. Their estimates are not far off 
Hubble's (1929) published value. 

To the present way of thinking the lengthy debate about the singularity in de Sitter's static solution, chronicled in North (1965), seems 
surprising, because de Sitter (1917) and Klein (1918) had presented de Sitter's solution as a sphere embedded in 4 plus 1 dimensional 
flat space, with no physical singularity. 



13 



is not affected by the error in the distance scale that affected Hubble's initial measurement of Ho). But the evidence 
now is that the mass density is about one quarter of what is predicted by this equation, as we will discuss. 

Einstein and de Sitter (1932) remark that the curvature term in Eq. jll] ) is "essentially determinable, and an 
increase in the precision of the data derived from observations will enable us in the future to fix its sign and determine 
its value." This is happening, 70 years later. The cosmological constant term is measurable in principle, too, and may 
now have been detected. But Einstein and de Sitter say only that the theory of an expanding universe with finite 
mean mass density "can be reached without the introduction of A. 

Further to this point, in the appendix of the second edition of his book, The Meaning of Relativity, Einstein (1945, 
p. 127) states that the "introduction of the 'cosmologic member"' — Einstein's terminology for A — "into the 
equations of gravity, though possible from the point of view of relativity, is to be rejected from the point of view 
of logical economy" , and that if "Hubble's expansion had been discovered at the time of the creation of the general 
theory of relativity, the cosmologic member would never have been introduced. It seems now so much less justified to 
introduce such a member into the field equations, since its introduction loses its sole original justification, — that of 
leading to a natural solution of the cosmologic problem." Einstein knew that without the cosmological constant the 
expansion time derived from Hubble's estimate of Hq is uncomfortably short compared to estimates of the ages of the 
stars, and opined that that might be a problem with the star ages. The big error, the value of Hq, was corrected by 
1960 (Sandage, 1958, 1962). 

Gamow (1970, p. 44) recalls that "when I was discussing cosmological problems with Einstein, he remarked that the 
introduction of the cosmological term was the biggest blunder he ever made in his life." This certainly is consistent 
with all of Einstein's written comments we have seen on the cosmological constant per se; we do not know whether 
Einstein was also referring to the missed chance to predict the evolution of the universe. 

B. The development of ideas 

1. Early indications of A 

In the classic book, The Classical Theory of Fields, Landau and Lifshitz (1951, p. 338) second Einstein's opinion 
of the cosmological constant A, stating there is "no basis whatsoever" for adjustment of the theory to include this 
term. The empirical side of cosmology is not much mentioned in this book, however (though there is a perceptive 
comment on the limited empirical support for the homogeneity assumption; p. 332). In the Supplementary Notes to 
the English translation of his book, Theory of Relativity, Pauli (1958, p. 220) also endorses Einstein's position. 

Discussions elsewhere in the literature on how one might find empirical constraints on the values of the cosmological 
parameters usually take account of A. The continued interest was at least in part driven by indications that A might 
be needed to reconcile theory and observations. Here are three examples. 

First, the expansion time is uncomfortably short if A = 0. Sandage's recalibration of the distance scale in the 1960s 
indicates Hq ~ 75 km s -1 Mpc -1 . If A = this says the time of expansion from densities too high for stars to have 
existed is < Hq 1 ~ 13 Gyr, maybe less than the ages of the oldest stars, then estimated to be greater than about 
15 Gyr. Sandage (1961a) points out that the problem is removed by adding a positive A. The present estimates 
reviewed below (Sec. IV. B. 3) are not far from these numbers, but still too uncertain for a significant case for A. 

Second, counts of quasars as a function of redshift show a peak at z ~ 2, as would be produced by the loitering epoch 
in Lemaitrc's A model (Pctrosian, Salpctcr, and Szekeres, 1967; Shklovsky, 1967; Kardashev, 1967). The peak is now 
well established, centered at z ~ 2.5 (Croom et al, 2001; Fan et al, 2001). It usually is interpreted as the evolution 
in the rate of violent activity in the nuclei of galaxies, though in the absence of a loitering epoch the indicated sharp 
variation in quasar activity with time is curious (but certainly could be a consequence of astrophysics that is not well 
understood). 

The third example is the redshift-magnitude relation. Sandage's (1961a) analysis indicates this is a promising 
method of distinguishing world models. The Gunn and Oke (1975) measurement of this relation for giant elliptical 
galaxies, with Tinsley's (1972) correction for evolution of the star population from assumed formation at high redshift, 
indicates curvature away from the linear relation in the direction that, as Gunn and Tinsley (1975) discuss, could 
only be produced by A (within general relativity theory). The new application of the redshift-magnitude test, to 
Type la supernovae (Sec. IV. B. 4), is not inconsistent with the Gunn-Oke measurement; we do not know whether this 
agreement of the measurements is significant, because Gunn and Oke were worried about galaxy evolution. 16 



Early measurements of the redshift-magnitude relation were meant in part to test the Steady State cosmology of Bondi and Gold (1948) 
and Hoyle (1948). Since the Steady State cosmology assumes spacetime is independent of time its line element has to have the form 
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2. The coincidences argument against A 

An argument against an observationally interesting value of A, from our distrust of accidental coincidences, has 
been in the air for decades, and became very influential in the early 1980s with the introduction of the inflation 
scenario for the very early universe. 

If the Einstein-de Sitter model in Eq. ([35|) were a good approximation at the present epoch, an observer measuring 
the mean mass density and Hubble's constant when the age of the universe was one tenth the present value, or at ten 
times the present age, would reach the same conclusion, that the Einstein-de Sitter model is a good approximation. 
That is, we would flourish at a time that is not special in the course of evolution of the universe. If on the other hand 
two or more of the terms in the expansion rate equation ( p"T| ) made substantial contributions to the present value of 
the expansion rate, it would mean we are present at a special epoch, because each term in Eq. ( pd] ) varies with the 
expansion factor in a different way. To put this in more detail, we imagine that the physics of the very early universe, 
when the relativistic cosmological model became a good approximation, set the values of the cosmological parameters. 
The initial values of the contributions to the expansion rate equation had to have been very different from each other, 
and had to have been exceedingly specially fixed, to make two of the fi^n's have comparable values. This would be 
a most remarkable and unlikely-looking coincidence. The multiple coincidences required for the near vanishing of d 
and a at a redshift not much larger than unity makes an even stronger case against Lemaitre's coasting model, by 
this line of argument. 

The earliest published comment we have found on this point is by Bondi (1960, p. 166), in the second edition of his 
book Cosmology. Bondi notes the "remarkable property" of the Einstein-de Sitter model: the dimensionless parameter 
we now call J7m is independent of the time at which it is computed (since it is unity). The coincidences argument 
follows and extends Bondi's comment. It is presented in McCrea (1971, p. 151). When Peebles was a postdoc, in the 
early 1960s, in R. H. Dicke's gravity research group, the coincidences argument was discussed, but published much 
later (Dicke, 1970, p. 62; Dicke and Peebles, 1979). We do not know its provenance in Dicke's group, whether from 
Bondi, McCrea, Dicke, or someone else. We would not be surprised to learn others had similar thoughts. 

The coincidences argument is sensible but not a proof, of course. The discovery of the 3 K thermal cosmic microwave 
background radiation gave us a term in the expansion rate equation that is down from the dominant one by four 
orders of magnitude, not such a large factor by astronomical standards. This might be counted as a first step away 
from the argument. The evidence from the dynamics of galaxies that JImo is less than unity is another step (Peebles, 
1984, p. 442; 1986). And yet another is the development of the evidence that the A and dark matter terms differ by 
only a factor of three (Eq. ||). This last is the most curious, but the community has come to accept it, for the most 
part. The precedent makes Lemaitre's coasting model more socially acceptable. 

A socially acceptable value of A cannot be such as to make life impossible, of course. 17 But perhaps the most 
productive interpretation of the coincidences argument is that it demands a search for a more fundamental underlying 
model. This is discussed further in Sec. III.E and the Appendix. 



3. Vacuum energy and A 

Another tradition to consider is the relation between A and the vacuum or dark energy density. In one approach 
to the motivation for the Einstein field equation, taken by McVittie (1956) and others, A appears as a constant of 
integration (of the expression for local conservation of energy and momentum). McVittie (1956, p. 35) emphasizes 
that, as a constant of integration, A "cannot be assigned any particular value on a priori grounds." Interesting 
variants of this line of thought are still under discussion (Weinberg, 1989; Unruh, 1989, and references therein). 

The notion of A as a constant of integration may be related to the issue of the zero point of energy. In laboratory 
physics one measures and computes energy differences. But the net energy matters for gravity physics, and one can 
imagine A represents the difference between the true energy density and the sum one arrives at by laboratory physics. 
Eddington (1939) and Lemaitre (1934, 1949) make this point. 

Bronstein (1933) 18 carries the idea further, allowing for transfer of energy between ordinary matter and that 



of the de Sitter solution with Qko = and the expansion parameter in Eq. (|27|). The measured curvature of the redshift- magnitude 
relation is in the direction predicted by the Steady State cosmology. But this cosmology fails other tests discussed in Sec. IV. B. 
If A were negative and the magnitude too large there would not be enough time for the emergence of life like us. If A were positive 
and too large the universe would expand too rapidly to allow galaxy formation. Our existence, which requires something resembling the 
Milky Way galaxy to contain and recycle heavy elements, thus provides an upper bound on the value of A. Such anthropic considerations 
are discussed by Weinberg (1987, 2001), Vilcnkin (2001), and references therein. 

Kragh (1996, p. 36) describes Bronstein's motivation and history. We discuss this model in more detail in Sec. III.E, and comment on 
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represented by A. In our notation, Bronstein expresses this picture by generalizing Eq. (g) to 

PA = -p-3-(p+p), (36) 
a 

where p and p are the energy density and pressure of ordinary matter and radiation. Bronstein goes on to propose 
a violation of local energy conservation, a thought that no longer seems interesting. North (1965, p. 81) finds 
Eddington's (1939) interpretation of the zero of energy also somewhat hard to defend. But for our purpose the 
important point is that the idea of A as a form of energy has been in the air, in at least some circles, for many years. 

The zero-point energy of fields contributes to the dark energy density. To make physical sense the sum over the 
zero-point mode energies must be cut off at a short distance or a high frequency up to which the model under 
consideration is valid. The integral of the zero-point energy (fc/2) of normal modes (of wavenumber k) of a massless 
real bosonic scalar field ($), up to the wavenumber cutoff k c , gives the vacuum energy density quantum-mechanical 
expectation value 19 

f k " 4irk 2 dk k k\ 
/>* = / "TOTo = ( 37 ) 
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Nernst (1916) seems to have been the first to write down this equation, in connection with the idea that the zero- 
point energy of the electromagnetic field fills the vacuum, as a light aether, that could have physically significant 
properties. 20 This was before Heisenberg and Schrodinger: Nernst's hypothesis is that each degree of freedom, which 
classical statistical mechanics assigns energy kT/2, has "Nullpunktsenergie" hu/2. This would mean the ground state 
energy of a one-dimensional harmonic oscillator is his, twice the correct value. Nernst's expression for the energy 
density in the electromagnetic field thus differs from Eq. (|37]) by a factor of two (after taking account of the two 
polarizations), which is wonderfully close. For a numerical example, Nernst noted that if the cutoff frequency were 
v = 10 20 Hz, or ~ 0.4 MeV, the energy density of the "Lichtather" (light aether) would be 10 23 erg cm -3 , or about 
100 g cm~ 3 . 

By the end of the 1920s Nernst's hypothesis was replaced with the demonstration that in quantum mechanics the 
zero-point energy of the vacuum is as real as any other. W. Pauli, in unpublished work in the 1920s, 21 repeated 
Nernst's calculation, with the correct factor of 2, taking k c to correspond to the classical electron radius. Pauli knew 
the value of pa is quite unacceptable: the radius of the static Einstein universe with this value of p\ "would not 
even reach to the moon" (Rugh and Zinkernagel, 2000, p. 5). 22 The modern version of this "physicists' cosmological 



why decay of dark energy into ordinary matter or radiation would be hard to reconcile with the thermal spectrum of the 3 K cosmic 
microwave background radiation. Decay into the dark sector may be interesting. 

Eq. (p7j), which usually figures in discussions of the vacuum energy puzzle, gives a helpful indication of the situation: the zero-point 
energy of each mode is real and the sum is large. The physics is seriously incomplete, however. The elimination of spatial momenta 
with magnitudes k > k c only makes sense if there is a preferred reference frame in which k c is defined. Magueijo and Smolin (2002) 
mention a related issue: In which reference frame is the Planck momentum of a virtual particle at the threshold for new phenomena? In 
both cases one may implicitly choose the rest frame for the large-scale distribution of matter and radiation. It seems strange to think 
the microphysics cares about large-scale structure, but maybe it happens in a sea of interacting fields. The cutoff in Eq. (jjTj) might 
be applied at fixed comoving wavenumber k c oc a(i) _1 , or at a fixed physical value of k c . The first prescription can be described by 
an action written as a sum of terms &^/% + k 2 &i / (2a(t) 2 ) for the allowed modes. The zero-point energy of each mode scales with the 

expansion of the universe as a(t) -1 , and the sum over modes scales as p$ oc a(t)~ 4 , consistent with k c oc a(t)~ 1 . In the limit of exact 
spatial homogeneity, an equivalent approach uses the spatial average of the standard expression for the field stress-energy tensor. Indeed, 
DeWitt (1975) and Akhmedov (2002) show that the vacuum expectation value of the stress-energy tensor, expressed as an integral cut 
off at k = k c , and computed in the preferred coordinate frame, is diagonal with space part pj> = p$/3, for the massless field we are 
considering. That is, in this prescription the vacuum zero point energy acts like a homogeneous sea of radiation. This defines a preferred 
frame of motion, where the stress-energy tensor is diagonal, which is not unexpected because we need a preferred frame to define k c . It 
is unacceptable as a model for the properties of dark energy, of course. For example, if the dark energy density were normalized to the 
value now under discussion, and varied as pj\ oc a(t)~ 4 , it would quite mess up the standard model for the origin of the light elements. 
We get a more acceptable model for the behavior of pa from the second prescription, with the cutoff at a fixed physical momentum. If 
we also want to satisfy local energy conservation we must take the pressure to be p,j = — p<j. This does not contradict the derivation of 
p§ in the first prescription, because the second situation cannot be described by an action: the pressure must be stipulated, not derived. 
What is worse, the known fields at laboratory momenta certainly do not allow this stipulation; they are well described by analogs of 
the action in the first prescription. This quite unsatisfactory situation illustrates how far we are from a theory of the vacuum energy. 
A helpful discussion of Nernst's ideas on cosmology is in Kragh (1996, pp. 151-7). 

This is discussed in Enz and Thellung (1960), Enz (1974), Rugh and Zinkernagel (2000, pp. 4-5), and Straumann (2002). 
In an unpublished letter in 1930, G. Gamow considered the gravitational consequences of the Dirac sea (Dolgov, 1989, p. 230). We 
thank A. Dolgov for helpful correspondence on this point. 
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constant problem" is even more acute, because a natural value for k c is thought to be much larger than what Nernst 
or Pauli used. 23 

While there was occasional discussion of this issue in the middle of the 20 th century (as in the quote from N. Bohr 
in Rugh and Zinkcrnagel, 2000, p. 5), the modern era begins with the paper by Zel'dovich (1967) that convinced the 
community to consider the possible connection between the vacuum energy density of quantum physics and Einstein's 
cosmological constant. 24 

If the physics of the vacuum looks the same to any inertial observer its contribution to the stress-energy tensor is 



the same as Einstein's cosmological constant (Eq. 19 1). Lemaitre (1934) notes this: "in order that absolute motion, 
i.e., motion relative to the vacuum, may not be detected, we must associate a pressure p = —pc 2 to the energy density 
pc 2 of vacuum". Gliner (1965) goes further, presenting the relation between the metric tensor and the stress-energy 
tensor of a vacuum that appears the same to any inertial observer. But it was Zel'dovich (1968) who presented the 
argument clearly enough and at the right time to catch the attention of the community. 

With the development of the concept of broken symmetry in the now standard model for particle physics came 
the idea that the expansion and cooling of the universe is accompanied by a sequence of first-order phase transitions 
accompanying the symmetry breaking. Each first-order transition has a latent heat that appears as a contribution 
to an effective time-dependent A(t) or dark energy density. 25 The decrease in value of the dark energy density at 
each phase transition is much larger than an acceptable present value (within relativistic cosmology); the natural 
presumption is that the dark energy is negligible now. This final condition seems bizarre, but the picture led to the 
very influential concept of inflation. We discussed the basic elements in connection with Eq. (p7j); we turn now to 
some implications. 



C. Inflation 



1. The scenario 



The deep issue inflation addresses is the origin of the large-scale homogeneity of the observable universe. In a 
relativistic model with positive pressure we can see distant galaxies that have not been in causal contact with each 
other since the singular start of expansion (Sec. II. C, Eq. |26|); they are said to be outside each other's particle horizon. 
Why do apparently causally unconnected parts of space look so similar? 26 Sato (1981a, 1981b), Kazanas (1980), and 
Guth (1981) make the key point: if the early universe were dominated by the energy density of a relatively flat real 
scalar field (inflaton) potential V(Q) that acts like A, the particle horizon could spread beyond the universe we can 
see. This would allow for the possibility that microphysics during inflation could smooth inhomogeneities sufficiently 
to provide an explanation of the observed large-scale homogeneity. (We are unaware of a definitive demonstration of 
this idea, however.) 

In the inflation scenario the field <J> rolls down its potential until eventually V(&) steepens enough to terminate 
inflation. Energy in the scalar field is supposed to decay to matter and radiation, heralding the usual Big Bang 
expansion of the universe. With the modifications of Guth's (1981) scenario by Linde (1982) and Albrecht and 



In terms of an energy scale ea denned by pa = e A , the Planck energy G -1 / 2 is about 30 orders of magnitude larger than the "observed" 
value of ea- This is, of course, an extreme case, since a lot of the theories of interest break down well below the Planck scale. Furthermore, 
in addition to other contributions, one is allowed to add a counterterm to Eq. ( [37] ) to predict any value of pa- With reference to this 
point, it is interesting to note that while Pauli did not publish his computation of pa, he remarks in his famous 1933 Handbuch der 
Physik review on quantum mechanics that it is more consistent to "exclude a zero-point energy for each degree of freedom as this energy, 
evidently from experience, does not interact with the gravitational field" (Rugh and Zinkernagel, 2000, p. 5). Pauli was fully aware 
that one must take account of zero-point energies in the binding energies of molecular structure, for example (and we expect he was 
aware that what contributes to the energy contributes to the gravitational mass). He chose to drop the section with the above comment 
from the second (1958) edition of the review (Pauli, 1980, pp. iv-v). In a globally supersymmetric field theory there are equal numbers 
of bosonic and fermionic degrees of freedom, and the net zero-point vacuum energy density p A vanishes (Iliopoulos and Zumino, 1974; 
Zumino, 1975). However, supcrsymmetry is not a symmetry of low energy physics, or even at the electroweak unification scale. It must 
be broken at low energies, and the proper setting for a discussion of the zero-point pa in this case is locally supersymmetric supergravity. 
Weinberg (1989, p. 6) notes "it is very hard to see how any property of supergravity or superstring theory could make the effective 
cosmological constant sufficiently small". Witten (2001) and Ellwanger (2002) review more recent developments on this issue in the 
superstring/M theory/branes scenario. 

For subsequent more detailed discussions of this issue, see Zel'dovich (1981), Weinberg (1989), Carroll, Press, and Turner (1992), Sahni 
and Starobinsky (2000), Carroll (2001), and Rugh and Zinkernagel (2000). 

Early references to this point are Linde (1974), Dreitlein (1974), Kirzhnitz and Linde (1974), Veltman (1975), Bludman and Ruderman 
(1977), Canuto and Lee (1977), and Kolb and Wolfram (1980). 

Early discussions of this question are reviewed by Rindler (1956); more recent examples are Misner (1969), Dicke and Peebles (1979), 
and Zee (1980). 
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Steinhardt (1982), the community quickly accepted this promising and elegant way to understand the origin of our 
homogeneous expanding universe. 27 

In Guth's (1981) picture the inflaton kinetic energy density is subdominant during inflation, <f> 2 <C V(&), so from 
Eqs. ( |bo|) the pressure p$ is very close to the negative of the mass density p$, and the expansion of the universe 
approximates the de Sitter solution, a oc exp(H\t) (Eq. p7|). 

For our comments on the spectrum of mass density fluctuations produced by inflation and the properties of solutions 
of the dark energy models in Sec. III.E we will find it useful to have another scalar field model. Lucchin and Matarrese 
(1985a, 1985b) consider the potential 

.4 



v($) = ^ ex P [-^V^qG\ , (38) 

where q and A are parameters. 28 They show that the scale factor and the homogeneous part of the scalar field evolve 
in time as 

o(i) = a [l + Nt} 2/q , 

$(t) - \n[l + Nt], (39) 



where N = 2qy/irA/ y/G(6 - q). If q < 2 this model inflates. Halliwell (1987) and Ratra and Peebles (1988) show 
that the solution ( p9| ) of the homogeneous equation of motion has the attractor property 29 mentioned in connection 
with Eq. ([nj). This exponential potential is of historical interest: it provided the first clear illustration of an attractor 
solution. We return to this point in Sec. III.E. 

A signal achievement of inflation is that it offers a theory for the origin of the departures from homogeneity. 
Inflation tremendously stretches length scales, so cosmologically significant lengths now correspond to extremely 
short lengths during inflation. On these tiny length scales quantum mechanics governs: the wavelenghts of zero-point 
field fluctuations generated during inflation are stretched by the inflationary expansion, 30 and these fluctuations are 
converted to classical density fluctuations in the late time universe. 31 

The power spectrum of the fluctuations depends on the model for inflation. If the expansion rate during inflation 
is close to exponential (Eq. pTj]), the zero-point fluctuations are frozen into primeval mass density fluctuations with 
power spectrum 

P(k) — {\8(k,t)\ 2 ) = AkT 2 {k). (40) 

Here 8{k, t) is the Fourier transform at wavenumber k of the mass density contrast 8(x, t) = p(x, t)/ (p(t)) — 1, where 
p is the mass density and (p) the mean value. After inflation, but at very large redshifts, the spectrum in this model 
is P(k) oc on all interesting length scales. This means the curvature fluctuations produced by the mass fluctuations 
diverge only as logfc. The form P(k) oc k thus need not be cut off anywhere near observationally interesting lengths, 
and in this sense it is scale-invariant. 32 The transfer function T(k) accounts for the effects of radiation pressure and 
the dynamics of nonrelativistic matter on the evolution of 8(k, t), computed in linear perturbation theory, at redshifts 
z /$. 10 4 . The constant A is determined by details of the chosen inflation model we need not get into. 
The exponential potential model in Eq. (^) produces the power spectrum 33 

P{k) = Ak n T 2 (k), n= (2-3q)/(2-q). (41) 

When n ^ l(g ^ 0) the power spectrum is said to be tilted. This offers a parameter n to be adjusted to fit the 
observations of large-scale structure, though as we will discuss the simple scale-invariant case n = 1 is close to the 
best fit to the observations. 



Aspects of the present state of the subject are reviewed by Guth (1997), Brandenberger (2001), and Lazarides (2002). 
Similar exponential potentials appear in some higher-dimensional Kaluza-Klein models. For an early discussion see Shan and Wetterich 
(1985). 

Ratra (1989, 1992a) shows that spatial inhomogeneities do not destroy this property, that is, for q < 2 the spatially inhomogeneous 
scalar field perturbation has no growing mode. 

The strong curvature of spacetime during inflation makes the vacuum state quite different from that of Minkowski spacetime (Ratra 
1985). This is somewhat analogous to how the Casimir metal plates modify the usual Minkowski spacetime vacuum state. 

For the development of these ideas see Hawking (1982), Starobinsky (1982), Guth and Pi (1982), Bardeen, Steinhardt, and Turner 
(1983), and Fischler, Ratra, and Susskind (1985). 

The virtues of a spectrum that is scale-invariant in this sense were noted before inflation, by Harrison (1970), Peebles and Yu (1970), 
and Zel'dovich (1972). 

This is discussed by Abbott and Wise (1984), Lucchin and Matarrese (1985a, 1985b), and Ratra (1989, 1992a). 
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The mass fluctuations in these inflation models are said to be adiabatic, because they are what you get by adiabat- 
ically compressing or decompressing parts of an exactly homogeneous universe. This means the initial conditions for 
the mass distribution are described by one function of position, 8(x,t). This function is a realization of a spatially 
stationary random Gaussian process, because it is frozen out of almost free quantum field fluctuations. Thus the sin- 
gle function of position is statistically prescribed by its power spectrum, as in Eqs. ( pio| ) and (pfl|). More complicated 
models for inflation produce density fluctuations that are not Gaussian, or do not have simple power law spectra, or 
have parts that break adiabaticity, as gravitational waves (Rubakov, Sazhin, and Veryaskin, 1982) or magnetic fields 
(Turner and Widrow, 1988; Ratra, 1992b) or new hypothetical fields. All these extra features may be invoked to fit 
the observations, if needed. It may be significant that none seem to be needed to fit the main cosmological structure 
constraints we have now. 



2. Inflation in a low density universe 

We do need an adjustment from the simplest case — an Einstein-de Sitter cosmology — to account for the mea- 
surements of the mean mass density. In the two models that lead to Eqs. ( fl(i| ) and pi] ) the enormous expansion 
factor during inflation suppresses the curvature of space sections, making Qko negligibly small. If A = 0, this fits the 
Einstein-de Sitter model (Eq. |35f] ), which in the absence of data clearly is the elegant choice. But the high mass 
density in this model was already seriously challenged by the data available in 1983, on the low streaming flow of the 
nearby galaxies toward the nearest known large mass concentration, in the Virgo cluster of galaxies, and the small 
relative velocities of galaxies outside the rich clusters of galaxies. 34 A striking and long familiar example of the latter 
is that the galaxies immediately outside the Local Group of galaxies, at distances of a few megaparsecs, are moving 
away from us in a good approximation to Hubble's homogeneous flow, despite the very clumpy distribution of galaxies 
on this scale. 35 The options (within general relativity) are that the mass density is low, so its clumpy distribution has 
little gravitational effect, or the mass density is high and the mass is more smoothly distributed than the galaxies. 
We comment on the first option here, and the second in connection with the cold dark matter model for structure 
formation in Sec. III.D. 

Under the first option we have two choices: introduce a cosmological constant, or space curvature, or maybe even 
both. In the conventional inflation picture space curvature is unacceptable, but there is another line of thought that 
leads to a universe with open space sections. Gott's (1982) scenario commences with a large energy density in an 
inflaton at the top of its potential. This behaves as Einstein's cosmological constant and produces a near de Sitter 
universe expanding asooc exp(H\t), with sufficient inflation to allow for a microphysical explanation of the large-scale 
homogeneity of the observed universe. As the inflaton gradually rolls down the potential it reaches a point where there 
is a small bump in the potential. The inflaton tunnels through this bump by nucleating a bubble. Symmetry forces 
the interior of the bubble to have open spatial sections (Coleman and De Luccia, 1980), and the continuing presence 
of a non-zero V(<f>) inside the bubble acts like A, resulting in an open inflating universe. The potential is supposed to 
steepen, bringing the second limited epoch of inflation to an end before space curvature has been completely redshifted 
away. The region inside the open bubble at the end of inflation is a radiation-dominated Friedmann-Lemaitre open 
model, with < Ukq < 1 (Eq. This can fit the dynamical evidence for low f^Mo with A = 0. 36 

The decision on which scenario, spatially-flat or open, is elegant, if either, depends ultimately on which Nature 
has chosen, if either. 37 But it is natural to make judgments in advance of the evidence. Since the early 1980s there 



This is discussed in Davis and Peebles (1983a, 1983b) and Peebles (1986). Relative velocities of galaxies in rich clusters are large, but 
the masses in clusters are known to add up to a modest mean mass density. Thus most of the Einstein-de Sitter mass would have to be 
outside the dense parts of the clusters, where the relative velocities are small. 

The situation a half century ago is illustrated by the compilation of galaxy redshifts in Humason, Mayall, and Sandage (1956). In this 
sample of 806 galaxies, 14 have negative redshifts (after correction for the rotation of the Milky Way galaxy and for the motion of the 
Milky Way toward the other large galaxy in the Local Group, the Andromeda Nebula), indicating motion toward us. Nine are members 
of the Local Group, at distances ^ 1 Mpc. Four are in the direction of the Virgo cluster, at redshift ~ 1200 km s and distance 
~ 20 Mpc. Subsequent measurements indicate two of these four really have negative redshifts, and plausibly are members of the Virgo 
cluster on the tail of the distribution of peculiar velocities of the cluster members. (Astronomers use the term peculiar velocity to denote 
the deviation from the uniform Hubble expansion velocity.) The last of the 14, NGC 3077, is in the M 81 group of galaxies at 3 Mpc 
distance. It is now known to have a small positive redshift. 

Gott's scenario is resurrected by Ratra and Peebles (1994, 1995). See Bucher and Turok (1995), Yamamoto, Sasaki, and Tanaka (1995), 
and Gott (1997), for further discussions of this model. In this case spatial curvature provides a second cosmologically-relevant length 
scale (in addition to that set by the Hubble radius H^ 1 ), so there is no natural preference for a power law power spectrum (Ratra, 
1994; Ratra and Peebles, 1995). 

At present, high energy physics considerations do not provide a compelling specific inflation model, but there are strong indications 
that inflation happens in a broad range of models, so it might not be unreasonable to think that future advances in high energy physics 
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have been occasional explorations of the open case, but the community generally has favored the flat case, f2xo = 0, 
without or, more recently, with a cosmological constant, and indeed the evidence now is that space sections are close 
to flat. The earlier preference for the Einstein-de Sitter case with JIko = and ^ao = led to considerable interest 
in the picture of biased galaxy formation in the cold dark matter model, as we now describe. 



D. The cold dark matter model 

Some of the present cosmological tests were understood in the 1930s; others are based on new ideas about structure 
formation. A decade ago a half dozen models for structure formation were under discussion 38 ; now the known viable 
models have been winnowed to one class: cold dark matter (CDM) and variants. We comment on the present state 
of tests of the CDM model in Sec. IV. A. 2, and in connection with the cosmological tests in Sec. IV. B. 

The CDM model assumes the mass of the universe now is dominated by dark matter that is nonbaryonic and acts 
like a gas of massive, weakly interacting particles with negligibly small primeval velocity dispersion. Structure is 
supposed to have formed by the gravitational growth of primeval departures from homogeneity that are adiabatic, 
scale- invariant, and Gaussian. The early discussions also assume an Einstein-de Sitter universe. These features all 
are naturally implemented in simple models for inflation, and the CDM model may have been inspired in part by 
the developing ideas of inflation. But the motivation in writing down this model was to find a simple way to show 
that the observed present-day mass fluctuations can agree with the growing evidence that the anisotropy of the 3 K 
thermal cosmic microwave background radiation is very small (Peebles, 1982). The first steps toward turning this 
picture into a model for structure formation were taken by Blumenthal et al. (1984). 

In the decade commencing about 1985 the standard cosmology for many active in research in this subject was 
the Einstein-de Sitter model, and for good reason: it eliminates the coincidences problem, it avoids the curiosity of 
nonzero dark energy, and it fits the condition from conventional inflation that space sections have zero curvature. But 
unease about the astronomical problems with the high mass density of the Einstein-de Sitter model led to occasional 
discussions of a low density universe with or without a cosmological constant, and the CDM model played an important 
role in these considerations, as we now discuss. 

When the CDM model was introduced it was known that the observations disfavor the high mass density of the 
Einstein-de Sitter model, unless the mass is more smoothly distributed than the visible matter (Sec. III.C). The key 
papers showing that this wanted biased distribution of visible galaxies relative to the distribution of all of the mass 
can follow in a natural way in the CDM theory are Kaiser (1984) and Davis et al. (1985). In short, where the mass 
density is high enough to lead to the gravitational assembly of a large galaxy the mass density tends to be high nearby, 
favoring the formation of neighboring large galaxies 

The biasing concept is important and certainly had to be explored. But in 1985 there was little empirical evidence 
for the effect and there were significant arguments against it, mainly the empty state of the voids between the 
concentrations of large galaxies. 39 In the biasing picture the voids contain most of the mass of an Einstein-de Sitter 
universe, but few of the galaxies, because galaxy formation there is supposed to have been suppressed. But it is 
hard to see how galaxy formation could be entirely extinguished: the CDM model would be expected to predict a 
void population of irregular galaxies, that show signs of a difficult youth. Many irregular galaxies are observed, but 
they avoid the voids. The straightforward reading of the observations thus is that the voids are empty, and that the 
dynamics of the motions of the visible galaxies therefore say Omo is well below unity, and that the mass is not more 
smoothly distributed than the visible galaxies. 

In a low density open universe, with = and positive ^ko, the growth of mass clustering is suppressed at 
z ^ fl^i —l. Thus to agree with the observed low redshift mass distribution density fluctuations at high redshift must 
be larger in the open model than in the Einstein-de Sitter case. This makes it harder to understand the small 3 K 
cosmic microwave background anisotropy. In a low density spatially-flat universe with f^Ko = and a cosmological 
constant, the transition from matter-dominated expansion to A-dominated expansion is more recent than the transition 



could give us a compelling and observationally successful model of inflation, that will determine whether it is flat or open. 
A scorecard is given in Peebles and Silk (1990). Structure formation models that assume all matter is baryonic, and those that augment 
baryons with hot dark matter such as low mass neutrinos, were already seriously challenged a decade ago. Vittorio and Silk (1985) show 
that the Uson and Wilkinson (1984) bound on the small-scale anisotropy of the 3 K cosmic microwave background temperature rules 
out a baryon-dominated universe with adiabatic initial conditions. This is because the dissipation of the baryon density fluctuations by 
radiation drag as the primeval plasma combines to neutral hydrogen (at redshift z ~ 1000) unacceptably suppresses structure formation 
on the scale of galaxies. Cold dark matter avoids this problem by eliminating radiation drag. This is one of the reasons attention turned 
to the hypothetical nonbaryonic cold dark matter. There has not been a thorough search for more baroque initial conditions that might 
save the baryonic dark matter model, however. 

The issue is presented in Peebles (1986); the data and history of ideas are reviewed in Peebles (2001). 
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from matter-dominated expansion to space-curvature-dominated expansion in an open universe with A = , as one 
sees from Eq. ([ll]). This makes density fluctuations grow almost as much as in the Einstein-de Sitter model, thus 
allowing smaller peculiar velocities in the flat-A case, a big help in understanding the observations. 40 

An argument for low JImoj with or without A, developed out of the characteristic length scale for structure in the 
CDM model. In the Friedmann-Lemaitre cosmology the mass distribution is gravitationally unstable. This simple 
statement has a profound implication: the early universe has to have been very close to homogeneous, and the growing 
departures from homogeneity at high redshift are well described by linear perturbation theory. The linear density 
fluctuations may be decomposed into Fourier components (or generalizations for open or closed space sections). At 
high enough redshift the wavelength of a mode is much longer than the time-dependent Hubble length H^ 1 , and 
gravitational instability makes the mode amplitude grow. Adiabatic fluctuations remain adiabatic, because different 
regions behave as if they were parts of different homogeneous universes. When the Hubble length becomes comparable 
to the mode proper wavelength, the baryons and radiation, strongly coupled by Thomson scattering at high redshift, 
oscillate as an acoustic wave and the mode amplitude for the cold dark matter stops growing. 41 The mass densities 
in dark matter and radiation are equal at redshift z cq — 2.4 x 10 4 i7Mo^ 2 ; thereafter the dark matter mass density 
dominates and the fluctuations in its distribution start to grow again. The suppressed growth of density fluctuations 
within the Hubble length at z > z Qq produces a bend in the power spectrum of the dark mass distribution, from 
P(k) oc k at long wavelengths, if scale-invariant (Eq. |EJ), to P(k) cx fc~ 3 at short wavelengths. 42 This means that 
at small scales (large k) the contribution to the variance of the mass density per logarithmic interval of wavelength 
is constant, and at small k the contribution to the variance of the Newtonian gravitational potential per logarithmic 
interval of k is constant. 

The wavelength at the break in the spectrum is set by the Hubble length at equal radiation and matter mass 
densities, t cq cx z~ q 2 . This characteristic break scale grows by the factor z cq to Abroak ~ z cq t cq at the present epoch. 
The numerical value is (Peebles, 1980a, Eq. [92.47]) 

Abrcak ~ 50n M lh~ 2 Mpc. (42) 

If A is close to constant, or A = and space sections are curved, it does not appreciably affect the expansion rate at 
redshift z cq , so this characteristic length is the same in flat and open cosmological models. In an Einstein-de Sitter 
model, with f^o = 1, the predicted length scale at the break in the power spectrum of CDM mass fluctuations is 
uncomfortably small relative to structures such as are observed in clusters of clusters of galaxies (superclusters), and 
relative to the measured galaxy two-point correlation function. That is, more power is observed on scales ~ 100 Mpc 
than is predicted in the CDM model with Qmo — 1- Since Abrcak scales as fij^o' a remedy is to go to a universe with 
small f^MOi either with A = and open space sections or f^xo = and a nonzero cosmological constant. The latter 
case is now known as ACDM. 43 



E. Dark energy 

The idea that the universe contains close to homogeneous dark energy that approximates a time- variable cosmolog- 
ical "constant" arose in particle physics, through the discussion of phase transitions in the early universe and through 
the search for a dynamical cancellation of the vacuum energy density; in cosmology, through the discussions of how 
to reconcile a cosmologically flat universe with the small mass density indicated by galaxy peculiar velocities; and on 
both sides by the thought that A might be very small now because it has been rolling toward zero for a very long 
time. 44 



The demonstration that the suppression of peculiar velocities is a lot stronger than the suppression of the growth of structure is in 
Peebles (1984) and Lahav et al. (1991). 

At high redshift the dark matter mass density is less than that of the radiation. The radiation thus fixes the expansion rate, which is 
too rapid for the self-gravity of the dark matter to have any effect on its distribution. Early discussions of this effect are in Guyot and 
Zel'dovich (1970) and Meszaros (1972) 

That is, the transfer function in Eq. (Ud) goes from a constant at small k to T 2 (fc) oc k~ 4 at large k. 

The scaling of A broak with Qmo was frequently noted in the 1980s. The earliest discussions we have seen of the significance for large- 
scale structure are in Silk and Vittorio (1987) and Efstathiou, Sutherland, and Maddox (1990), who consider a spatially-flat universe, 
Blumenthal, Dekel, and Primack (1988), who consider the open case, and Holtzman (1989), who considers both. For the development 
of tests of the open case see Lyth and Stewart (1990), Ratra and Peebles (1994), Kamionkowski et al. (1994a), and Gorski et al. 
(1995). Pioneering steps in the analysis of the anisotropy of the 3 K cosmic microwave background temperature in the ACDM model 
include Kofman and Starobinsky (1985) and Gorski, Silk, and Vittorio (1992). We review developments after the COBE detection of 
the anisotropy (Smoot et al., 1992) in Sees. IV.B.ll and 12. 

This last idea is similar in spirit to Dirac's (1937, 1938) attempt to explain the large dimensionless numbers of physics. He noted 
that the gravitational force between two protons is much smaller than the electromagnetic force, and that that might be because the 
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The idea that the dark energy is decaying by emission of matter or radiation is now strongly constrained by the 
condition that the decay energy must not significantly disturb the spectrum of the 3 K cosmic microwave background 
radiation. But the history of the idea is interesting, and decay to dark matter still a possibility, so we comment on 
both here. The picture of dark energy in the form of defects in cosmic fields has not received much attention in recent 
years, in part because the computations are difficult, but might yet prove to be productive. Much discussed nowadays 
is dark energy in a slowly varying scalar field. The idea is reviewed at some length here and in even more detail in 
the Appendix. We begin with another much discussed approach: prescribe the dark energy by parameters in numbers 
that seem fit for the quality of the measurements. 



1. The XCDM parametrization 

In the XCDM parametrization the dark energy interacts only with itself and gravity, the dark energy density 
Px(t) > is approximated as a function of world time alone, and the pressure is written as 

Px = w x px, (43) 

an expression that has come to be known as the cosmic equation of state. 45 Then the local energy conservation 
equation (||) is 

|=-3-px(l + »x). (44) 
at a 

If wx is constant the dark energy density scales with the expansion factor as 

px«^ 3(1+Mx) . (45) 

If wx < — 1/3 the dark energy makes a positive contribution to a/a (Eq. ||). If wx — — 1/3 the dark energy has 
no effect on a, and the energy density varies as « 1/a 2 , the same as the space curvature term in a 2 /a 2 (Eq. PI). 
That is, the expansion time histories are the same in an open model with no dark energy and in a spatially-flat model 
with wx = —1/3, although the spacetime geometries differ. 46 If wx < —1 the dark energy density is increasing. 47 

Equation (Eq) with constant wx has the great advantage of simplicity. An appropriate generalization for the more 
precise measurements to come might be guided by the idea that the dark energy density is close to homogeneous, 
spatial variations rearranging themselves at close to the speed of light, as in the scalar field models discussed below. 
Then for most of the cosmological tests we have an adequate general description of the dark energy if we let wx be a 
free function of time. In scalar field pictures wx is derived from the field model; it can be a complicated function 
of time even when the potential is a simple function of the scalar field. 



gravitational constant G is decreasing in inverse proportion to the world time. This is the earliest discussion we know of what has come 
to be called the hierarchy problem, that is, the search for a mechanism that might be responsible for the large ratio between a possibly 
more fundamental high energy scale, for example, that of grand unification or the Planck scale (where quantum gravitational effects 
become significant) and a lower possibly less fundamental energy scale, for example that of electroweak unification (see, for example, 
Georgi, Quinn, and Weinberg, 1974). The hierarchy problem in particle physics may be rephrased as a search for a mechanism to 
prevent the light electroweak symmetry breaking Higgs scalar field mass from being large because of a quadratically divergent quantum 
mechanical correction (see, for example, Susskind, 1979). In this sense it is similar in spirit to the physicists' cosmological constant 
problem of Sec. III.B. 

Other parametrizations of dark energy are discussed by Hu (1998) and Bucher and Spergel (1999). The name, XCDM, for the case 
wx < in Eq. (fiil), was introduced by Turner and White (1997). There is a long history in cosmology of applications of such an 
equation of state, and the related evolution of p/±; examples are Canuto et al. (1977), Lau (1985), Huang (1985), Fry (1985), Hiscock 
(1986), Ozer and Taha (1986), and Olson and Jordan (1987). See Ratra and Peebles (1988) for references to other early work on a 
time-variable A and Overduin and Cooperstock (1998) and Sahni and Starobinsky (2000) for reviews. More recent discussions of this 
and related models may be found in John and Joseph (2000), Zimdahl et al. (2001), Dalai et al. (2001), Gudmundsson and Bjornsson 
(2002), Bean and Melchiorri (2002), Mak, Belinchon, and Harko (2002), and Kujat et al. (2002), through which other recent work may 
be traced. 

As discussed in Sec. IV, it appears difficult to reconcile the case tox = — 1/3 with the Type la supernova apparent magnitude data 
(Garnavich et al., 1998; Perlmutter et al., 1999a). 

This is quite a step from the thought that the dark energy density is small because it has been rolling to zero for a long time, but the 
case has found a context (Caldwell, 2002; Maor et al., 2002). Such models were first discussed in the context of inflation (e.g., Lucchin 
and Matarrese, 1985b), where it was shown that the u>x < — 1 component could be modeled as a scalar field with a negative kinetic 
energy density (Peebles, 1989a). 

The availability of a free function greatly complicates the search for tests as opposed to curve fitting! This is clearly illustrated by Maor 
et al. (2002). For more examples see Perlmutter, Turner, and White (1999b) and Efstathiou (1999). 
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The analysis of the large-scale anisotropy of the 3 K cosmic microwave background radiation requires a prescription 
for how the spatial distribution of the dark energy is gravitationally related to the inhomogeneous distribution of 
other matter and radiation (Caldwell et al., 1998). In XCDM this requires at least one more parameter, an effective 
speed of sound, with c 2 x > (for stability, as discussed in Sec. II. B), in addition to wx- 



2. Decay by emission of matter or radiation 

Bronstein (1933) introduced the idea that the dark energy density p\ is decaying by the emission of matter or 
radiation. The continuing discussions of this and the associated idea of decaying dark matter (Sciama, 2001, and 
references therein) are testimony to the appeal. Considerations in the decay of dark energy include the effect on the 
formation of light elements at z ~ 10 10 , the contribution to the 7-ray or optical extragalactic background radiation, 
and the perturbation to the spectrum of the 3 K cosmic microwave background radiation. 

The effect on the 3 K cosmic microwave background was of particular interest a decade ago, as a possible explanation 
of indications of a significant departure from a Planck spectrum. Precision measurements now show the spectrum is 
very close to thermal. The measurements and their interpretation are discussed by Fixsen et al. (1996). They show 
that the allowed addition to the 3 K cosmic microwave background energy density pr is limited to just Sps./ pn £ 10~ 4 
since redshift z ~ 10 5 , when the interaction between matter and radiation was last strong enough for thermal 
relaxation. The bound on 5pn/ pn is not inconsistent with what the galaxies are thought to produce, but it is well 
below an observationally interesting dark energy density. 

Dark energy could decay by emission of dark matter, cold or hot, without disturbing the spectrum of the 3 K cosmic 
microwave background radiation. For example, let us suppose the dark energy equation of state is wx = — 1, and 
hypothetical microphysics causes the dark energy density to decay as pa a ~ n by the production of nonrelativistic 
dark matter. Then Bronstein's Eq. (^6|) says the dark matter density varies with time as 

PM(t) = ^3 + 33^/>a(*), (46) 

where A is a constant and < n < 3. In the late time limit the dark matter density is a fixed fraction of the dark 
energy. But for the standard interpretation of the measured anisotropy of the 3 K background we would have to 
suppose the first term on the right hand side of Eq. (Efj|) is not much smaller than the second, so the coincidences 
issue discussed in Sec. III.B.2 is not much relieved. It does help relieve the problem with the small present value of 
PA (to be discussed in connection with Eq. Q ) . 

We are not aware of any work on this decaying dark energy picture. Attention instead has turned to the idea that 
the dark energy density evolves without emission, as illustrated in Eq. ( ftSTj and the two classes of physical models to 
be discussed next. 



3. Cosmic field defects 

The physics and cosmology of topological defects produced at phase transitions in the early universe are reviewed by 
Vilenkin and Shcllard (1994). An example of dark energy is a tangled web of cosmic string, with fixed mass per unit 
length, which self-intersects without reconnection. In Vilenkin's (1984) analysis 50 the mean mass density in strings 
scales as String oc (ta(t)) . When ordinary matter is the dominant contribution to a 2 /a 2 , the ratio of mass densities 
is Pstring/p °c t 1 / 3 . Thus at late times the string mass dominates. In this limit, ^string oc a" 2 , Wx = — l/3 for the 
XCDM parametrization of Eq. (^5j), and the universe expands as a cx t. Davis (1987) and Kamionkowski and Toumbas 
(1996) propose the same behavior for a texture model. One can also imagine domain walls fill space densely enough 



These considerations generally are phenomenological: the evolution of the dark energy density, and its related coupling to matter or 
radiation, is assigned rather than derived from an action principle. Recent discussions include Pollock (1980), Kazanas (1980), Freese et 
al. (1987), Gasperini (1987), Sato, Terasawa, and Yokoyama (1989), Bartlett and Silk (1990), Overduin, Wesson, and Bowyer (1993), 
Matyjasek (1995), and Birkel and Sarkar (1997). 

The string flops at speeds comparable to light, making the coherence length comparable to the expansion time t. Suppose a string 
randomly walks across a region of physical size a(t)R in N steps, where aR ~ N 1 ' 2 t. The total length of this string within the region R 
is I ~ Nt. Thus the mean mass density of the string scales with time as p s tring <x l/a? oc (ta(t))~ 1 . One randomly walking string does 
not fill space, but we can imagine many randomly placed strings produce a nearly smooth mass distribution. Spergel and Pen (1997) 
compute the 3 K cosmic microwave background radiation anisotropy in a related model, where the string network is fixed in comoving 
coordinates so the mean mass density scales as p a tring oc a -2 . 
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not to be dangerous. If the domain walls are fixed in comoving coordinates the domain wall energy density scales as 
px oc a" 1 (Zel'dovich, Kobzarev, and Okun, 1974; Battye, Bucher, and Spergel, 1999). The corresponding equation 
of state parameter is wx — —2/3, which is thought to be easier to reconcile with the supernova measurements than 
wx = —1/3 (Garnavich et at, 1998; Perhnutter et ah, 1999a). The cosmological tests of defects models for the dark 
energy have not been very thoroughly explored, at least in part because an accurate treatment of the behavior of the 
dark energy is difficult (as seen, for example, in Spergel and Pen, 1997; Friedland, Muruyama, and Perelstein, 2002), 
but this class of models is worth bearing in mind. 



4. Dark energy scalar field 

At the time of writing the popular picture for dark energy is a classical scalar field with a self-interaction potential 
V(<&) that is shallow enough that the field energy density decreases with the expansion of the universe more slowly 
than the energy density in matter. This idea grew in part out of the inflation scenario, in part from ideas from particle 
physics. Early examples are Weiss (1987) and Wetterich (1988). 51 The former considers a quadratic potential with 
an ultralight effective mass, an idea that reappears in Frieman et al. (1995). The latter considers the time variation 
of the dark energy density in the case of the Lucchin and Matarrese (1985a) exponential self-interaction potential 
(Eq. @). 52 

In the exponential potential model the scalar field energy density varies with time in constant proportion to the 
dominant energy density. The evidence is that radiation dominates at redshifts in the range 10 3 ^ z ^ 10 10 , from the 
success of the standard model for light element formation, and matter dominates at 1 ;$ z ^ 10 3 , from the success of 
the standard model for the gravitational growth of structure. This would leave the dark energy subdominant today, 
contrary to what is wanted. This led to the proposal of the inverse power-law potential in Eq. ([3l]) for a single real 
scalar field. 53 

We do not want the hypothetical field $ to couple too strongly to baryonic matter and fields, because that would 
produce a "fifth force" that is not observed. 5 ' 55 Within quantum field theory the inverse power-law scalar field 
potential makes the model non-renormalizable and thus pathological. But the model is meant to describe what might 
emerge out of a more fundamental quantum theory, which maybe also resolves the physicists' cosmological constant 
problem (Sec. III.B), as the effective classical description of the dark energy. 56 The potential of this classical effective 



51 Other early examples include those cited in Ratra and Peebles (1988) as well as Endo and Fukui (1977), Fujii (1982), Dolgov (1983), 
Nilles (1985), Zee (1985), Wilczek (1985), Bertolami (1986), Ford (1987), Singh and Padmanabhan (1988), and Barr and Hochberg 
(1988). 

52 For recent discussions of this model see Ferreira and Joyce (1998), Ott (2001), Hwang and Noh (2001), and references therein. 

53 In what follows we focus on this model, which was proposed by Peebles and Ratra (1988). The model assumes a conventionally 
normalized scalar field kinetic energy and spatial gradient term in the action, and it assumes the scalar field is coupled only to itself 
and gravity. The model is then completely characterized by the form of the potential (in addition to all the other usual cosmological 
parameters, including initial conditions). Models based on other forms for V(&), with a more general kinetic energy and spatial gradient 
term, or with more general couplings to gravity and other fields, are discussed in the Appendix. 

54 The current value of the mass associated with spatial inhomogeneities in the field is m^(to) ~ Ho ~ 10 — 33 eV, as one would expect from 
the dimensions. More explicitly, one arrives at this mass by writing the field as <t>(i, x) = (<E>) (i) + <f>(t , x) and Taylor expanding the scalar 
field potential energy density V($) about the homogeneous mean background (<t>) to quadratic order in the spatially inhomogeneous part 
<j>, to get m 2 = V" (($)). Within the context of the inverse power-law model, the tiny value of the mass follows from the requirements 
that V varies slowly with the field value and that the current value of V be observationally acceptable. The difference between the roles 
of m,£ and the constant m q in the quadratic potential model V = m 2 <I> 2 /2 is worth noting. The mass m q has an assigned and arguably 
fine-tuned value. The effective mass mj, ~ H belonging to V oc $ — a is a derived quantity, that evolves as the universe expands. The 
small value of m^(to) explains why the scalar field energy cannot be concentrated with the non-relativistic mass in galaxies and clusters 
of galaxies. Because of the tiny mass a scalar field would mediate a new long-range fifth force if it were not weakly coupled to ordinary 
matter. Weak coupling also ensures that the contributions to coupling constants (such as the gravitational constant) from the exchange 
of dark energy bosons are small, so such coupling constants are not significantly time variable in this model. See, for example, Carroll 
(1998), Chiba (1999), Horvat (1999), Amendola (2000), Bartolo and Pietroni (2000), and Fujii (2000) for recent discussions of this and 
related issues. 

55 Coupling between dark energy and dark matter is not constrained by conventional fifth force measurements. An example is discussed 
by Amendola and Tocchini-Valentini (2002). Perhaps the first consideration is that the fifth- force interaction between neighboring dark 
matter halos must not be so strong as to shift regular galaxies of stars away from the centers of their dark matter halos. 

56 Of course, the zero-point energy of the quantum-mechanical fluctuations around the mean field value contributes to the physicists' 
cosmological constant problem, and renormalization of the potential could destroy the attractor solution (however, see Doran and Jackel, 
2002) and could generate couplings between the scalar field and other fields leading to an observationally inconsistent "fifth force". The 
problems within quantum field theory with the idea that the energy of a classical scalar field is the dark energy, or drives inflation, are 
further discussed in the Appendix. The best we can hope is that the effective classical model is a useful approximation to what actually 
is happening, which might lead us to a more satisfactory theory. 
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field is chosen ad hoc, to fit the scenario. But one can adduce analogs within supergravity, superstring/M, and brane 
theory, as reviewed in the Appendix. 

The solution for the mass fraction in dark energy in the inverse power-law potential model (in Eq. |33| when p$ <C p, 
and the numerical solution at lower redshifts) is not unique, but it behaves as what has come to be termed an attractor 
or tracker: it is the asymptotic solution for a broad range of initial conditions. 57 The solution also has the property 
that p$ is decreasing, but less rapidly than the mass densities in matter and radiation. This may help alleviate two 
troubling aspects of the cosmological constant. The coincidences issue is discussed in Sec. III.B. The other is the 
characteristic energy scale set by the value of A, 

ca(*o) = Pa(£o) 1/4 = 0.003(1 - ft M0 ) 1/4 /i 1/2 eV, (47) 

when f2j^ an d ^ko may be neglected. In the limit of constant dark energy density, cosmology seems to indicate new 
physics at an energy scale more typical of chemistry. If p\ is rolling toward zero the energy scale might look more 
reasonable, as follows (Peebles and Ratra, 1988; Steinhardt et al., 1999; Brax et al., 2000). 

Suppose that as conventional inflation ends the scalar field potential switches over to the inverse power-law form in 
Eq. (|3~i|). Let the energy scale at the end of inflation be e(ti) = p(ii) 1 / 4 , where p(t\) is the energy density in matter 
and radiation at the end of inflation, and let CA(ii) be the energy scale of the dark energy at the end of inflation. 
Since the present value e\(to) of the dark energy scale (Eq. J47fl) is comparable to the present energy scale belonging 
to the matter, we have from Eq. ( |33| ) 

e. U W 2/(« + 2 ) 



*M*«v{M) • (48) 

For parameters of common inflation models, e(ti) ~ 10 13 GeV, and e\(to)/e(tj) ~ 1CP 25 . If, say, a = 6, then 

e A (ti) ~ 10~ 6 e(O - 10 7 GeV. (49) 

As this example illustrates, one can arrange the scalar field model so it has a characteristic energy scale that exceeds 
the energy ~ 10 3 GeV below which physics is thought to be well understood: in this model cosmology does not force 
upon us the idea that there is as yet undiscovered physics at the very small energy in Eq. (|47j). Of course, where the 
factor ~ 10~ 6 in Eq. ( (49| ) comes from still is an open question, but, as discussed in the Appendix, perhaps easier to 
resolve than the origin of the factor ~ 1CP 25 in the constant A case. 

When we can describe the dynamics of the departure from a spatially homogeneous field in linear perturbation 
theory, a scalar field model generally is characterized by the time-dependent values of wx (Eq. [Q) and the speed of 
sound c s x (e.g., Ratra, 1991; Caldwell et al., 1998). In the inverse power-law potential model the relation between 
the power-law index a and the equation of state parameter in the matter-dominated epoch is independent of time 
(Ratra and Quillen, 1992), 

w x = 2 —. (50) 

a + 2 

When the dark energy density starts to make an appreciable contribution to the expansion rate the parameter wx 
starts to evolve. The use of a constant value of wx to characterize the inverse power-law potential model thus can 
be misleading. For example, Podariu and Ratra (2000, Fig. 2) show that, when applied to the Type la supernova 
measurements, the XCDM parametrization in Eq. (|50|) leads to a significantly tighter apparent upper limit on wx, at 
fixed f^MOj than is warranted by the results of a computation of the evolution of the dark energy density in this scalar 
field model. Caldwell et al. (1998) deal with the relation between scalar field models and the XCDM parametrization 
by fixing wx, as a constant or some function of rcdshift, deducing the scalar field potential V($) that produces this 
wx, and then computing the gravitational response of $ to the large-scale mass distribution. 



A recent discussion is in Brax and Martin (2000). Brax, Martin, and Riazuelo (2000) present a thorough analysis of the evolution 
of spatial inhomogeneities in the inverse power-law scalar field potential model and confirm that these inhomogencities do not destroy 
the homogeneous attractor solution. For other recent discussions of attractor solutions in a variety of contexts see Liddle and Scherrer 
(1999), Uzan (1999), de Ritis et al. (2000), Holden and Wands (2000), Baccigalupi, Matarrese, and Perrotta (2000), and Huey and 
Tavakol (2002). 
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IV. THE COSMOLOGICAL TESTS 

Our intention is to supplement recent discussions of parameter determinations within the standard relativistic 
cosmology 58 with a broader consideration of the issues summarized in two questions: what is the purpose of the 
cosmological tests, and how well is the purpose addressed by recent advances and work in progress? 

The short answer to the first question used to be that we seek to check the underlying physical theory, general 
relativity, applied on the time and length scales of cosmology; the model for the stress-energy tensor in Einstein's field 
equation, suitably averaged over the rich small-scale structure we cannot describe in any detail; and the boundary 
condition, that the universe we can observe is close to homogeneous and isotropic on the scale of the Hubble length. 
Recent advances make use of the CDM prescription for the stress-energy tensor and the boundary condition, so we 
must add the elements of the CDM model to the physics to be checked. 

The short answer to the second question is that we now have searching checks of the standard cosmology, which 
the model passes. But we believe it takes nothing away from the remarkable advances of the tests, and the exemplary 
care in the measurements, to note that there is a lot of room for systematic errors. As we discussed in Sec. LA, the 
empirical basis for the standard model for cosmology is not nearly as substantial as is the empirical basis for the 
standard model for particle physics: in cosmology it is not yet a matter of measuring parameters in a well-established 
physical theory. 

We comment on the two main pieces of physics, general relativity and the CDM model, in Sees. IV. A. 1 and IV. A. 2. 
In Sec. IV. B we discuss the state of 13 cosmological tests, proceeding roughly in order of increasing model dependence. 
We conclude that there is a well-established scientific case for the physical significance of the matter density parameter, 
and for the result of the measurements, 0.15 ^ Quo ^ 0.4 (in the sense of a two standard deviation range). Our 
reasoning is summarized in Sec. IV. C, along with an explanation of why we are not so sure about the detection of A 
or dark energy. 

A. The theories 

1. General relativity 

Some early discussions of the cosmological tests, as in Robertson (1955) and Bondi (1960), make the point that 
observationally important elements of a spatially homogeneous cosmology follow by symmetry, independent of general 
relativity. This means some empirical successes of the cosmology are not tests of relativity. The point was important 
in the 1950s, because the Steady State theory was a viable alternative to the Friedmann-Lemaitre cosmology, and 
because the experimental tests of relativity were quite limited. 

The tests of general relativity are much better now, but cosmology still is a considerable extrapolation. The length 
scales characteristic of the precision tests of general relativity in the Solar System and binary pulsar are ^ 10 13 cm. 
An important scale for cosmology is the Hubble length, Hq 1 ~ 5000 Mpc ~ 10 28 cm, fifteen orders of magnitude 
larger. An extrapolation of fifteen orders of magnitude in energy from that achieved at the largest accelerators, 
~ 10 12 eV, brings us to the very different world of the Planck energy. Why is the community not concerned about 
an extrapolation of similar size in the opposite direction? One reason is that the known open issues of physics have 
to do with small length scales; there is no credible reason to think general relativity may fail on large scales. This 
is comforting, to be sure, but, as indicated in footnote 7, not the same as a demonstration that we really know the 
physics of cosmology. Another reason is that if the physics of cosmology were very different from general relativity it 
surely would have already been manifest in serious problems with the cosmological tests. This also is encouraging, 
but we have to consider details, as follows. 

One sobering detail is that in the standard cosmology the two dominant contributions to the stress-energy tensor — 
dark energy and dark matter — are hypothetical, introduced to make the theories fit the observations (Eq. ||). This 
need not mean there is anything wrong with general relativity — we have no reason to expect Nature to have made 
all matter readily observable other than by its gravity — but it is a cautionary example of the challenges. Milgrom's 
(1983) modified Newtonian dynamics (MOND) replaces the dark matter hypothesis with a hypothetical modification of 
the gravitational force law. MOND gives remarkably successful fits to observed motions within galaxies, without dark 
matter (de Blok et ai, 2001). So why should we believe there really is cosmologically significant mass in nonbaryonic 
dark matter? Unless we are lucky enough to get a laboratory detection, the demonstration must be through the 
tests of the relativistic cosmology (and any other viable cosmological models that may come along, perhaps including 



See Bahcall et at (1999), Schindler (2001), Sarkar (2002), Freedman (2002), Plionis (2002), and references therein. 
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an extension of MOND). This indirect chain of evidence for dark matter is becoming tight. A new example — the 
prospect for a test of the inverse square law for gravity on the length scales of cosmology — is striking enough for 
special mention here. 59 

Consider the equation of motion 60 of a freely moving test particle with nonrelativistic peculiar velocity v in a 
universe with expansion factor a(t), 

dv a _ . 1 . . 

— + -v = g = —V<p. (51) 
at a a 

The particle always is moving toward receding observers, which produces the second term in the left-most expression. 
The peculiar gravitational acceleration g relative to the homogeneous background model is computed from the Poisson 
equation for the gravitational potential <p, 

\7 2 ip = 4TrGa 2 [p(x,t)-(p)}. (52) 

The mean mass density (p) is subtracted because g is computed relative to the homogeneous model. The equation 
of mass conservation expressed in terms of the density contrast S = p/(p) — 1 of the mass distribution modeled as a 
continuous pressureless fluid is 

c)$ 1 

— + -V- (1 + 5)v = Q. (53) 
at a 

In linear perturbation theory in v and S these equations give 

W + 2 ~a% = i7TG{p)6 ' ^ t] = f & D ®- (54) 

Here D(t) is the growing solution to the first equation. 61 The velocity field belonging to the solution D(t) is the 
inhomogeneous solution to Eq. ( [53] ) in linear perturbation theory, 



* = n-Av) *v, f - ^ 6 o- (55) 



Ijha f y-x f/: ( .^ o ,,„ 
47r J \y-x\ 3 

The factor f — d log D j d log a depends on the cosmological model; the second equation is a good approximation 
if A = or space curvature vanishes. 62 One sees from Eq. ([35]) that the peculiar velocity is proportional to the 
gravitational acceleration, as one would expect in linear theory. 

The key point of Eq. ( pi} ) for the present purpose is that the evolution of the density contrast 5 at a given position 
is not affected by the value of S anywhere else. This is a consequence of the inverse square law. The mass fluctuation 
in a chosen volume element produces a peculiar gravitational acceleration Sg that produces a peculiar velocity field 
Sv oc g that has zero divergence and so the mass inside the volume element does not effect the exterior. 

For a "toy" model of the effect of a failure of the inverse square law, suppose we adjust the expression for the 
peculiar gravitational acceleration produced by a given mass distribution to 

g = a 3 R [ d 3 yS(y)-^^-Q(a\y- x\), (56) 



\y- 

where R is some function of world time only. In standard gravity physics Q(w) = w~ 2 . We have no basis in 
fundamental physics for any other function of w. Although Milgrom's (1983) MOND provides a motivation, Eq. j5^ ) 
is not meant to be an extension of MOND to large-scale flows. It is an ad hoc model that illustrates an important 
property of the inverse square law. 



Binetruy and Silk (2001) and Uzan and Bernardeau (2001) pioneered this probe of the inverse square law. Related probes, based on 
the relativistic dynamics of gravitational lensing and the anisotropy of the 3 K thermal background, are discussed by these authors and 
White and Kochanek (2001). 

These relations are discussed in many books on cosmology, including Peebles (1980a). 

The general solution is a sum of the growing and decaying solutions, but because the universe has expanded by a large factor since 
nongravitational forces were last important on large scales we can ignore the decaying part. 

This is illustrated in Fig. 13.14 in Peebles (1993). An analytic expression for spherical symmetry is derived by Lightman and Schechter 
(1990). 
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We noted that in linear theory v oc g. Thus we find that the divergence of Eq. ( pq ) , with the mass conservation 
equation ( |53|) in linear perturbation theory, gives 

^ + 2 ^ =S{k,t)S n , S(k,t) = (4xRk/a)J w 2 dwQ{w)j x {kw/a), (57) 

where <%(i) is the Fourier transform of the mass density contrast 5(x,t) and ji is a spherical Bessel function. The 
inverse square law, Q = w -2 , makes the factor S independent of the wavenumber k. This means all Fourier amplitudes 
grow by the same factor in linear perturbation theory (when the growing mode dominates), so the functional form of 
S(x, t) is preserved and the amplitude grows, as Eq. (]54j) says. When Q is some other function, the phases of the rJg 
are preserved but the functional form of the power spectrum \8^\ 2 evolves. For example, if Q oc w n ~ 2 with —2 < n < 1 
(so the integral in Eq. j57| does not diverge) Eq. ( p7| ) is 

<9 2 <% _ a dSz 



dt 2 



where [/ is some function of world time. 

If n > density fluctuations grow faster on larger scales. If Q(w) follows Newtonian gravity on the scale of galaxies 
and bends to n > on larger scales it reduces the mean mass density needed to account for the measured large-scale 
galaxy flows, and maybe reduces the need for dark matter. But there are testable consequences: the apparent value 
of S^mo would vary with the length scale of the measurement, and the form of the power spectrum of the present 
mass distribution would not agree with the form at redshift z ~ 1000 when it produced the observed angular power 
spectrum of the 3 K cosmic microwave background. Thus we are very interested in the evidence of consistency of 
these tests (as discussed in Sec. IV.B.13). 



2. The cold dark matter model for structure formation 



Important cosmological tests assume the CDM model for structure formation (Sec. III.C), so we must consider 
tests of the model. The model has proved to be a useful basis for analyses of the physics of formation of galaxies and 
clusters of galaxies (e.g., Kay et ai, 2002; Colberg et ai, 2000; and references therein). There are issues to consider, 
however; Sellwood and Kosowsky (2001) give a useful survey of the situation. We remark on recent developments and 
what seem to us to be critical issues. 

Numerical simulations of the dark mass distribution in the CDM model predict that massive halos have many low 
mass satellites, perhaps significantly more than the number observed around the Milky Way galaxy (Klypin et ai, 
1999; Moore et ai, 1999a). The issue is of great interest but not yet a critical test, because it is difficult to predict 
the nature of star formation in a low mass dark halo: what does a dark halo look like when star formation or the 
neutral gas content makes it visible? For recent discussions see Tully et al. (2002) and Stoehr et al. (2002). 

The nature of the dark mass distribution within galaxies is a critical issue, because we know where to look for a 
distinctive CDM feature: a cusp-like central mass distribution, the density varying with radius r as p oc r~ a with 
a 1. The power law is not unexpected, because there is nothing in the CDM model to fix an astronomically 
interesting value for a core radius. 63 A measure of the mass distribution in disk galaxies is the rotation curve: the 
circular velocity as a function of radius for matter supported by rotation. In some low surface brightness galaxies the 
observed rotation curves are close to solid body, v c oc r, near the center, consistent with a near homogeneous core, 
and inconsistent with the cusp- like CDM mass distribution. 64 

The circular velocity produced by the mass distribution p oc r _1 is not very different from solid body, or from the 
observations, and the difference might be erased by gravitational rearrangement of the dark mass by the fluctuations 
in the distribution of baryonic mass driven by star formation, winds, or supernovae. This is too complicated to assess 
by current numerical simulations. But we do have a phenomenological hint: central solid body rotation is most clearly 
seen in the disk-like galaxies with the lowest surface brightnesses, the objects in which the baryon mass seems least 
likely to have had a significant gravitational effect on the dark mass. This challenge to the CDM model is pressing. 



Pioneering work on the theory of the central mass distribution in a dark mass halo is in Dubinski and Carlberg (1991). Moore (1994) 
and Flores and Primack (1994) are among the first to point out the apparent disagreement between theory and observation. 

The situation is reviewed by de Blok et al. (2001), and de Blok and Bosma (2002). The galaxy NGC 3109 is a helpful example because 
it is particularly close — just outside the Local Group — and so particularly well resolved. An optical image is in plate 39 in the Hubble 
Atlas of Galaxies (Sandage, 1961b). The radial velocity measurements across the face of the galaxy, in Figs. 1 and 2 in Blais-Ouellette, 
Amram, and Carignan (2001), are consistent with circular motion with v c oc r at r ,$ 2 kpc. 
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The challenge may be resolved in a warm dark matter model, where the particles are assigned a primeval velocity 
dispersion that suppresses the initial power spectrum of density fluctuations on small scales (Moore et al., 1999b; 
Sommer-Larsen and Dolgov, 2001; Bode, Ostriker, and Turok, 2001). But it seems to be difficult to reconcile the 
wanted suppression of small-scale power with the observation of small-scale clustering in the Lyman-a forest — the 
neutral hydrogen observed at z ~ 3 in the Lyman-a resonance absorption lines in quasar spectra (Narayanan et al., 
2000; Knebe et al, 2002). Spergel and Stcinhardt (2000) point out that the scattering cross section of self-interacting 
cold dark matter particles can be adjusted to suppress the cusp-like core. 65 Dave et al. (2001) demonstrate the effect 
in numerical simulations. But Miralda-Escude (2002) points out that the collisions would tend to make the velocity 
distribution isotropic, contrary to the evidence for ellipsoidal distributions of dark matter in clusters of galaxies. For 
recent surveys of the very active debate on these issues see Primack (2002) and Tasitsiomi (2002); for references to 
still other possible fixes see Dave et al. (2001). 

Another critical issue traces back to the biasing picture discussed in Sec. III.D. If Omo is well below unity there need 
not be significant mass in the voids defined by the large galaxies. But the biasing process still operates, and might 
be expected to cause dwarf or irregular galaxies to trespass into the voids outlined by the large regular galaxies. This 
seems to happen in CDM model simulations to a greater extent than is observed. Mathis and White (2002) discuss 
voids in ACDM simulations, but do not address the trespassing issue. The reader is invited to compare the relative 
distributions of big and little galaxies in the simulation in Fig. 1 of Mathis and White (2002) with the examples of 
observed distributions in Figs. 1 and 2 in Peebles (1989b) and in Figs. 1 to 3 in Peebles (2001). 

The community thought is that the trespassing issue need not be a problem for the CDM model: the low mass density 
in voids disfavors formation of galaxies from the debris left in these regions. But we have not seen an explanation of 
why local upward mass fluctuations, of the kind that produce normal galaxies in populated regions, and appear also 
in the predicted debris in CDM voids, fail to produce dwarf or irregular void galaxies. An easy explanation is that the 
voids contain no matter, having been gravitationally emptied by the growth of primeval non-Gaussian mass density 
fluctuations. The evidence in tests (10) and (11) in Sec. IV. B is that the initial conditions are close to Gaussian. But 
non-Gaussian initial conditions that reproduce the character of the galaxy distribution, including suppression of the 
trespassing effect, would satisfy test (10) by construction. 

We mention finally the related issues of when the large elliptical galaxies formed and when they acquired the central 
compact massive objects that are thought to be remnant quasar engines (Lynden-Bell, 1969). 

In the CDM model large elliptical galaxies form in substantial numbers at redshift z < 1. Many astronomers do 
not see this as a problem, because ellipticals do tend to contain relatively young star populations, and some elliptical 
galaxies have grown by recent mergers, as predicted in the CDM model. 66 But prominent merger events are rare, 
and the young stars seen in ellipticals generally seem to be a "frosting" (Trager et al., 2000) of recent star formation 
on a dominant old star population. The straightforward reading of the evidence assembled in Peebles (2002) is that 
most of the large ellipticals are present as assembled galaxies of stars at z = 2. 67 The ACDM model prediction is 
uncertain because it depends on the complex processes of star formation that are so difficult to model. The reading 
of the situation by Thomas and Kauffmann (1999) is that the predicted abundance of giant ellipticals at z — 2 is less 
than about one third of what it is now. Deciding whether the gap between theory and observation can be closed is 
not yet straightforward. 

A related issue is the significance of the observations of quasars at redshift z ~ 6. By conventional estimates 68 



In a power law halo with p <x r -7 , the velocity dispersion varies with radius as (v 2 ) ~ GM(< r)/r oc r 2— 7 . The particle scattering cross 
section must be adjusted to erase the effective temperature gradient, thus lowering the mass density at small radii, without promoting 
unacceptable core collapse. 

The classic merger example is also the nearest large elliptical galaxy, Ccntaurus A (NGC 5128). The elliptical image is crossed by 
a band of gas and dust that likely is the result of a merger with one of the spiral galaxies in the group around this elliptical. For a 
thorough review of what is known about this galaxy see Israel (1998). 

Papovich, Dickinson, and Ferguson (2002) find evidence that the comoving number density of all galaxies with star mass greater than 
1 X 10 10 M Q , where M Q is the mass of the Sun, is significantly less at redshift z > 1 than now. This is at least roughly in line with the 
distribution of star ages in the Milky Way spiral galaxy: the bulge stars are old, while the stars in the thin disk have a broad range of 
ages. Thus if this galaxy evolved from z = 2 without significant growth by mergers its star mass at z = 2 would be significantly less 
than the present value, which is about 5 X 1O 1O M0. Cimatti et al. (2002) show that the redshift distribution of faint galaxies selected 
at wavelenght A ~ 2 (im is not inconsistent with the picture that galaxy evolution at z < 2 is dominated by ongoing star formation 
rather than merging. 

The quasars discovered in the Sloan Digital Sky Survey are discussed by Fan et al. (2001). If the quasar radiation is not strongly 
beamed toward us, its luminosity translates to an Eddington mass (the mass at which the gravitational pull on unshielded plasma 
balances the radiation pressure) M bll ~ 1O 9,3 M0. In a present-day elliptical galaxy with this mass in the central compact object the 
line of sight velocity dispersion is <r ~ 350 km s _1 . This is close to the highest velocity dispersion observed in low redshift elliptical 
galaxies. For example, in the Faber et al. (1989) catalog of 500 ellipticals, 15 have 300 < a < 400 km s _1 , and none has a larger a. 
From the present-day relation between a and luminosity, an elliptical galaxy with a = 350 km s _1 has mass ~ 1O 12 3 M in stars. The 
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these quasars are powered by black holes with masses at the upper end of the range of masses of the central compact 
objects — let us call them black hole quasar remanants — in the largest present-day elliptical galaxies. Here are 
some options to consider. First, the high redshift quasars may be in the few large galaxies that have already formed 
at z ~ 6. Wyithe and Loeb (2002), following Efstathiou and Rees (1988), show that this fits the ACDM model if 
the quasars at z ~ 6 have black hole mass ~ 10 9 M Q in dark halos with mass <~ 10 12 M Q . In the ACDM picture 
these early galaxies would be considerably denser than normal galaxies; to be checked is whether they would be rare 
enough to be observationally acceptable. Second, the quasars at z <~ 6 may be in more modest star clusters that later 
grew by merging into giant ellipticals. To be established is whether this growth would preserve the remarkably tight 
correlation between the central black hole mass and the velocity dispersion of the stars 69 , and whether growth by 
merging would produce an acceptable upper bound on black hole masses at the present epoch. Third, large ellipticals 
might have grown by accretion around pre-existing black holes, without a lot of merging. This is explored by Danese 
et al. (2002). 

There does not seem to be a coherent pattern to the present list of challenges to the CDM model. The rotation 
curves of low surface brightness galaxies suggest we want to suppress the primeval density fluctuations on small 
scales, but the observations of what seem to be mature elliptical galaxies at high redshifts suggest we want to increase 
small-scale fluctuations, or maybe postulate non-Gaussian fluctuations that grow into the central engines for quasars 
at z ~ 6. We do not want these central engines to appear in low surface brightness galaxies, of course. 

It would not be at all surprising if the confusion of challenges proved to be at least in part due to the difficulty 
of comparing necessarily schematic analytic and numerical model analyses to the limited and indirect empirical 
constraints. But it is also easy to imagine that the CDM model has to be refined because the physics of the dark 
sector of matter and energy is more complicated than ACDM, and maybe even more complicated than any of the 
alternatives now under discussion. Perhaps some of the structure formation ideas people were considering a decade ago, 
which invoke good physics, also will prove to be significant factors in relieving the problems with structure formation. 
And the important point for our purpose is that we do not know how the relief might affect the cosmological tests. 



B. The tests 

The literature on the cosmological tests is enormous compared to what it was just a decade ago, and growing. Our 
references to this literature are much sparser than in Sec. Ill, on the principle that no matter how complete the list it 
will be out of date by the time this review is published. For the same reason, we do not attempt to present the best 
values of the cosmological parameters based on their joint fit to the full suite of present measurements. The situation 
will continue to evolve as the measurements improve, and the state of the art is best followed on astro-ph. We do 
take it to be our assignment to consider what the tests are testing, and to assess the directions the results seem to 
be leading us. The latter causes us to return many times to two results that seem secure because they are so well 
checked by independent lines of evidence, as follows. 

First, at the present state of the tests, optically selected galaxies are useful mass tracers. By that we mean the 
assumption that visible galaxies trace mass does not seriously degrade the accuracy of analyses of the observations. 
This will change as the measurements improve, of course, but the case is good enough now that we suspect the 
evidence will continue to be that optically selected galaxies are good indicators of where most of the mass is at the 
present epoch. Second, the mass density in matter is significantly less than the critical Einstein-de Sitter value. The 
case is compelling because it is supported by so many different lines of evidence (as summarized in Sec. IV. C). Each 
could be compromised by systematic error, to be sure, but it seems quite unlikely the evidence could be so consistent 
yet misleading. A judgement of the range of likely values of the mass density is more difficult. Our estimate, based 
on the measurements we most trust, is 

0.15 < Qmo £ 0.4, (59) 

and we would put the central value at 12mo — 0.25. The spread is meant in the sense of two standard deviations: we 
would be surprised to find f^Mo is outside this range. 



dark matter associated with this many baryons is Mdm ~ 10 Mq. This is a large mass to assemble at z ~ 6, but it helps that such 
objects are rare. The present number density of giant elliptical galaxies with c > 300 km s" 1 is about 10 -5 Mpc -3 , four orders of 
magnitude more than the comoving number density of quasars detected at z ~ 6. 

Ferrarese and Merritt (2000) and Gebhardt et al. (2000) show that the black hole mass correlates with the velocity dispersion of the 
stars in an elliptical galaxy and the velocity dispersion of the bulge stars in a spiral galaxy. This is not a direct gravitational effect: the 
black hole mass is less than 1% of the star mass in the bulge or the elliptical galaxy. 
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Several other policy decisions should be noted. The first is that we do not comment on tests that have been 
considered but not yet applied in a substantial campaign of measurements. A widely discussed example is the Alcock 
and Paczyhski (1979) comparison of the apparent depth and width of a system from its angular size and depth in 
redshift. 

In analyses of the tests of models for evolving dark energy density, simplicity recommends the XCDM parametriza- 
tion with a single constant parameter wx, as is demonstrated by the large number of recent papers on this approach. 
But the more complete physics recommends the scalar field model with an inverse power-law potential. This includes 
the response of the spatial distribution of the dark energy to the peculiar gravitational field. Thus our comments 
on variable dark energy density are more heavily weighted to the scalar field model than is the case in the recent 
literature. 

The gravitational deflection of light appears not only as a tool in cosmological tests, as gravitational lensing, but 
also as a source of systematic error. The gravitational deflections caused by mass concentrations magnify the image 
of a galaxy along a line of sight where the mass density is larger than the average, and reduce the solid angle of the 
image when the mass density along the line of sight is low. The observed energy flux density is proportional to the 
solid angle (because the surface brightness, erg cm~ 2 s _1 ster -1 Hz -1 , is conserved at fixed redshift). Selection can be 
biased either way, by the magnification effect or by obscuration by the dust that tends to accompany mass. 70 When 
the tests are more precise we will have to correct them for these biases, through models for the mass distribution 
(as in Premadi et at, 2001), and the measurements of the associated gravitational shear of the shapes of the galaxy 
images. But the biases seem to be small and will not be discussed here. 

And finally, as the cosmological tests improve a satisfactory application will require a joint fit of all of the parameters 
to all of the relevant measurements and constraints. Until recently it made sense to impose prior conditions, most 
famously the hope that if the universe is not well described by the Einstein-de Sitter model then surely it is the 
case either that A is negligibly small or else that space curvature may be neglected. We suspect the majority in the 
community still expect this is true, on the basis of the coincidences argument in Sec. II. B. 2, but it will be important 
to see what comes out of joint fits of both Omo and Oao, as well as all the other parameters, as is becoming the 
current practice. Our test-by-test discussion is useful for sorting out the physics and astronomy, we believe; it is not 
the prototype for the coming generations of precision application of the tests. 

Our remarks are ordered by our estimates of the model dependence. 



1. The thermal cosmic microwave background radiation 

We are in a sea of radiation with spectrum very close to Planck at T = 2.73 K, and isotropic to one part in 10 5 
(after correction for a dipole term that usually is interpreted as the result of our motion relative to the rest frame 
defined by the radiation). The thermal spectrum indicates thermal relaxation, for which the optical depth has to 
be large on the scale of the Hubble length Hq 1 . We know space now is close to transparent at the wavelengths of 
this radiation, because radio galaxies are observed at high redshift. Thus the universe has to have expanded from a 
state quite different from now, when it was hotter, denser, and optically thick. This is strong evidence our universe 
is evolving. 

This interpretation depends on, and checks, conventional local physics with a single metric description of spacetime. 
Under these assumptions the expansion of the universe preserves the thermal spectrum and cools the temperature 
as 72 

Toe (! + *). (60) 



This was recognized by Zel'dovich (1964), R. Feynman, in 1964, and S. Refsdal, in 1965. Feynman's comments in a colloquium are noted 
by Gunn (1967). Peebles attended Refsdal's lecture at the International Conference on General Relativity and Gravitation, London, 
July 1965; Refsdal (1970) mentions the lecture. 

The history of the discovery and measurement of this radiation, and its relation to the light element abundances in test (2), is presented 
in Peebles (1971, pp. 121-9 and 240-1), Wilkinson and Peebles (1990), and Alpher and Herman (2001). The precision spectrum 
measurements are summarized in Halpern, Gush, and Wishnow (1991) and Fixsen et al. (1996). 

To see this, recall the normal modes argument used to get Eq. (pp. The occupation number in a normal mode with wavelength A at 
temperature T is the Planck form A/* = [e^ c /' sTA — 1] — Adiabaticity says M is constant. Since the mode wavelength varies as A oc a(i), 
where a is the expansion factor in Eq. (kh, and Af is close to constant, the mode temperature varies as T oc l/a(t). Since the same 
temperature scaling applies to each mode, an initially thermal sea of radiation remains thermal in the absence of interactions. We do 
not know the provenance of this argument; it was familiar in Dicke's group when the 3 K cosmic microwave background radiation was 
discovered. 
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Bahcall and Wolf (1968) point out that one can test this temperature-redshift relation by measurements of the 
excitation temperatures of fine-structure absorption line systems in gas clouds along quasar lines of sight. The 
corrections for excitations by collisions and the local radiation field are subtle, however, and perhaps not yet fully 
sorted out (as discussed by Molaro et al., 2002, and references therein). 

The 3 K thermal cosmic background radiation is a centerpiece of modern cosmology, but its existence does not test 
general relativity. 



2. Light element abundances 



The best evidence that the expansion and cooling of the universe traces back to high redshift is the success of 
the standard model for the origin of deuterium and isotopes of helium and lithium, by reactions among radiation, 
leptons, and atomic nuclei as the universe expands and cools through temperature T ~ 1 MeV at redshift z ~ 10 10 . 
The free parameter in the standard model is the present baryon number density. The model assumes the baryons 
are uniformly distributed at high redshift, so this parameter with the known present radiation temperature fixes the 
baryon number density as a function of temperature and the temperature as a function of time. The latter follows 



from the expansion rate Eq. (11), which at the epoch of light element formation may be written as 

-V = !W, (6i) 

a J 3 

where the mass density p r counts radiation, which is now at T = 2.73 K, the associated neutrinos, and e pairs. The 
curvature and A terms are unimportant, unless the dark energy mass density varies quite rapidly. 

Independent analyses of the fit to the measured element abundances, corrected for synthesis and destruction in 
stars, by Buries, Nollett, and Turner (2001), and Cyburt, Fields, and Olive (2001), indicate 

0.018 < VtBoh 2 < 0.022, and 0.006 < fl B oh 2 < 0.017, (62) 

both at 95% confidence limits. Other analyses by Coc et al. (2002) and Thuan and Izotov (2002) result in ranges 
that lie between the two of Eq. ([32]). The difference in values may be a useful indication of remaining uncertainties; 
it is mostly a consequence of the choice of isotopes used to constrain f^Bo^ 2 - Buries et al. (2001) use the deuterium 
abundance, Cyburt et al. (2001) favor the helium and lithium measurements, and the other two groups use other 
combinations of abundances. Equation. (|6^ ) is consistent with the summary range, 0.0095 < ^Boh 2 < 0.023 at 95% 
confidence, of Fields and Sarkar (2002). 

The baryons observed at low redshift, in stars and gas, amount to (Fukugita, Hogan, and Peebles, 1998) 

^bo ~ 0.01. (63) 

It is plausible that the difference between Eqs. d62| ) and ( |63| ) is in cool plasma, with temperature T ~ 100 eV, in 
groups of galaxies. It is difficult to observationally constrain the idea that there is a good deal more cool plasma in 
the large voids between the concentrations of galaxies. A more indirect but eventually more precise constraint on 
fiBOj from the anisotropy of the 3 K thermal cosmic microwave background radiation, is discussed in test (11). 

It is easy to imagine complications, such as inhomogeneous entropy per baryon, or in the physics of neutrinos; 
examples may be traced back through Abazajian, Fuller, and Patel (2001) and Giovannini, Keihanen, and Kurki- 
Suonio (2002). It seems difficult to imagine that a more complicated theory would reproduce the successful predictions 
of the simple model, but Nature fools us on occasion. Thus before concluding that the theory of the pre-stellar light 
element abundances is known, apart from the addition of decimal places to the cross sections, it is best to wait and 
see what advances in the physics of baryogenesis and of neutrinos teach us. 

How is general relativity probed? The only part of the computation that depends specifically on this theory is the 
pressure term in the active gravitational mass density, in the expansion rate equation (|§|) . If we did not have general 
relativity, a simple Newtonian picture might have led us to write down a/a = —47rGp r /3 instead of Eq. (S). With 
p r ~ 1/a 4 , as appropriate since most of the mass is fully relativistic at the redshifts of light element production, this 
would predict the expansion time a/a is 2 1 / 2 times the standard expression (that from Eq. |6l| ). The larger expansion 
time would hold the neutron to proton number density ratio close to that at thermal equilibrium, n/p = e~®' kT , 
where Q is the difference between the neutron and proton masses, to lower temperature. It would also allow more time 
for free decay of the neutrons after thermal equilibrium is broken. Both effects decrease the final 4 He abundance. The 
factor 2 1 / 2 increase in expansion time would reduce the helium abundance by mass to Y ~ 0.20. This is significantly 
less than what is observed in objects with the lowest heavy element abundances, and so seems to be ruled out 
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(Steigman, 20 02). 73 That is, we have positive evidence for the relativistic expression for the active gravitational mass 
density at redshift z ~ 10 10 , a striking result. 



3. Expansion times 



The predicted time of expansion from the very early universe to redshift z is 



da f°° dz 



where -E(z) is defined in Eq. (ftl]). If A = the present age is to < Hq 1 . In the Einstein-de Sitter model the present 
age is to = 2/(3ifn). If the dark energy density is significant and evolving, we may write p\ = pAof( z )i where the 
function of redshift is normalized to /(0) = 1. Then E{z) generalizes to 

E{z) = (O M0 (l + zf + n m (l + zf + n K0 (l + + n AQ f(z)) 1/2 . (65) 

In the XCDM parametrization with constant wx (Eq. {f5|), f(z) = (1 + z) 3 ( 1+tux ). Olson and Jordan (1987) present 
the earliest discussion we have found of iToto in this picture (before it got the name). In scalar field models, f(z) 
generally must be evaluated numerically; examples are in Peebles and Ratra (1988). 

The relativistic correction to the active gravitational mass density (Eq. |j) is not important at the redshifts at 
which galaxies can be observed and the ages of their star populations estimated. At moderately high redshift, where 
the nonrelativistic matter term dominates, Eq. ( |64| ) is approximately 

t{z)^—^{l + z)-^. (66) 

That is, the ages of star populations at high redshift are an interesting probe of ^mo but they are not very sensitive 
to space curvature or to a near constant dark energy density. 74 

Recent analyses of the ages of old stars 75 indicate the expansion time is in the range 

11 Gyr < t < 17 Gyr, (67) 

at 95% confidence, with central value to — 13 Gyr. Following Krauss and Chaboyer (2001) these numbers add 0.8 Gyr 
to the star ages, under the assumption star formation commenced no earlier than z = 6 (Eq. |66f| ). A naive addition 
in quadrature to the uncertainty in Hq (Eq. ||) indicates the dimensionless age parameter is in the range 

0.72 < H t £ 1.17, (68) 

at 95% confidence, with central value Hoto ~ 0.89. The uncertainty here is dominated by that in to. In the spatially- 
flat ACDM model (Oko = 0), Eq. (|6|) translates to 0.15 ;$ Qmo ~ 0.8, with central value Qmo — 0.4. In the open 
model with J7ao = 0, the constraint is ^mo ~ 0-6 with the central value ^mo — 0.1. In the inverse power-law scalar 
field dark energy case (Sec. II. C) with power-law index a = 4, the constraint is 0.05 £ I^mo ~ 0.8. 

We should pause to admire the unification of the theory and measurements of stellar evolution in our galaxy which 
yield the estimate of to, and the measurements of the extragalactic distance scale, which yield Hq, in the product in 



Eq. (68) that agrees with the relativistic cosmology with dimensionless parameters in the range now under discussion. 
As we indicated in Sec. Ill, there is a long history of discussion of the expansion time as a constraint on cosmological 
models. The measurements now are tantalizingly close to a check of consistency with the values of Omo and f^Ao 
indicated by other cosmological tests. 



There is a long history of discussions of this probe of the expansion rate at the redshifts of light element production. The reduction of 
the helium abundance to Y ~ 0.2 if the expansion time is increased by the factor 2 1 / 2 is seen in Figs. 1 and 2 in Peebles (1966). Dicke 
(1968) introduced the constraint on evolution of the strength of the gravitational interaction; see Uzan (2002) for a recent review. The 
effect of the number of neutrino families on the expansion rate and hence the helium abundance is noted by Hoyle and Tayler (1964) 
and Shvartsman (1969). Steigman, Schramm, and Gunn (1977) discuss the importance of this effect as a test of cosmology and of the 
particle physics measures of the number of neutrino families. 

The predicted maximum age of star populations in galaxies at redshifts 2^1 does still depend on f2A.o and CkO; an d there is the 
advantage that the predicted maximum age is a lot shorter than today. This variant of the expansion time test is discussed by Nolan 
et al. (2001), Lima and Alcaniz (2001), and references therein. 

See Carretta et al. (2000), Krauss and Chaboyer (2001), Chaboyer and Krauss (2002), and references therein. 
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4. The redshift-angular size and redshift-magnitude relations 

An object at redshift z with physical length I perpendicular to the line of sight subtends angle 6 such that 

l = a(t)r(z)6 = a r(z)9/(l + z), (69) 

where clq = a(to). The angular size distance r(z) is the coordinate position of the object in the first line element in 
Eq. (|l5|), with the observer placed at the origin. The condition that light moves from source to observer on a radial 
null geodesic is 

r(z) dr f dt 

(70) 



VI + Kr 2 J a(t)' 



which gives 

dz 

V^ko V Jo E(z 



H a r(z) = sinh (0^0 / ^TT ) - ( 71 ) 

'K0 V Jo 



where E(z) is defined in Eqs. (yj) and (pq). 

In the Einstein-de Sitter model, the angular-size-redshift relation is 

a H i(i+ Z y/ 2 



2 VTTI-i 



(72) 



At z <C 1, ^ = HqI/z, consistent with the Hubble redshift-distance relation. At z 1 the image is magnified, 76 
oc I + z. 

The relation between the luminosity of a galaxy and the energy flux density received by an observer follows from 
Liouville's theorem: the observed energy flux i Vo per unit time, area, solid angle, and frequency satisfies 

= i v Jv e /(\ + z) 4 , (73) 

with i Ue the emitted energy flux (surface brightness) at the source and Sv e = 5iso(l + z) the bandwidth at the source 
at redshift z. The redshift factor (I + z) 4 appears for the same reason as in the 3 K cosmic microwave background 
radiation energy density. With Eq. ([39]) to fix the solid angle, Eq. (jT^) says the observed energy flux per unit area, 
time, and frequency from a galaxy at redshift z that has luminosity L Vc per frequency interval measured at the source 
is 

f "° = 47rogr(z)2(l + 2)' (74) 

In conventional local physics with a single metric theory the redshift-angular size (Eq. |6l| ) and redshift-magnitude 
(Eq. [j74j) relations are physically equivalent. 77 

The best present measurement of the redshift-magnitude relation uses supernovae of Type la. 78 The results are 
inconsistent with the Einstein-de Sitter model, at enough standard deviations to make it clear that unless there is 



The earliest discussion we know of the magnification effect is by Hoyle (1959). In the coordinate system in Eq. jl5j), with the observer 
at the origin, light rays from the object move to the observer along straight radial lines. An image at high redshift is magnified because 
the light detected by the observer is emitted when the proper distance to the object measured at fixed world time is small. Because 
the proper distance between the object and source is increasing faster than the speed of light, emitted light directed at the observer is 
initially moving away from the observer. 

For a review of measurements of the redshift-magnitude relation (and other cosmological tests) we recommend Sandage (1988). A recent 
application to the most luminous galaxies in clusters is in Aragon-Salamanca, Baugh, and Kauffmann (1998). The redshift-angular 
size relation is measured by Daly and Guerra (2001) for radio galaxies, Buchalter et al. (1998) for quasars, and Gurvits, Kellermann, 
and Frey (1999) for compact radio sources. Constraints on the cosmological parameters from the Gurvits et al. data are discussed 
by Vishwakarma (2001), Lima and Alcaniz (2002), Chen and Ratra (2003), and references therein, and constraints based on the radio 
galaxy data are discussed by Daly and Guerra (2001), Podariu et al. (2003), and references therein. 
These supernovae are characterized by the absence of hydrogen lines in the spectra; they are thought to be the result of explosive 
nuclear burning of white dwarf stars. Pskovskii (1977) and Phillips (1993) pioneered the reduction of the supernovae luminosities to a 
near universal standard candle. For recent discussions of their use as a cosmological test see Goobar and Perlmutter (1995), Reiss et al. 
(1998), Perlmutter et al. (1999a), Gott et al. (2001), and Leibundgut (2001). We recommend Leibundgut's (2001) cautionary discussion 
of astrophysical uncertainties: the unknown nature of the trigger for the nuclear burning, the possibility that the Phillips correction to 
a fiducial luminosity actually depends on redshift or environment within a galaxy, and possible obscuration by intergalactic dust. There 
are also issues of physics that may affect this test (and others): the strengths of the gravitational or electromagnetic interactions may 
vary with time, and photon-axion conversion may reduce the number of photons reaching us. All of this is under active study. 
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something quite substantially and unexpectedly wrong with the measurements the Einstein-de Sitter model is ruled 
out. The data require A > at two to three standard deviations, depending on the choice of data and method of 
analysis (Leibundgut, 2001; Gott et al., 2001). The spatially-flat case with fluo in the range of Eq. ( |59| ) is a good 
fit for constant A. The curre nt data do not provi de interesting constraints on the models for evolving dark energy 
density. Perlmutter et al. (http://snap.lbl.gov/) show that a tighter constraint, from supernovae observations to 
redshift z ~ 2, by the proposed SNAP satellite, is feasible and capable of giving a significant detection of A and maybe 
its evolution. 80 



5. Galaxy counts 

Counts of galaxies — or of other objects whose number density as a function of redshift may be modeled — probe 
the volume element (dV/dz)dzSQ defined by a solid angle SD, in the sky and a redshift interval dz. The volume is 
fixed by the angular size distance (Eq. p3|), which determines the area subtended by the solid angle, in combination 
with the redshift-time relation (Eq. |p4[| jTwhich fixes the radial distance belonging to the redshift interval. 

Sandage (1961a) and Brown and Tinsley (1974) showed that with the technology then available galaxy counts are 
not a very sensitive probe of the cosmological parameters. Loh and Spillar (1986) opened the modern exploration 
of the galaxy count-redshift relation at redshifts near unity, where the predicted counts are quite different in models 
with and without a cosmological constant (as illustrated in Figure 13.8 in Peebles, 1993). 

The interpretation of galaxy counts requires an understanding of the evolution of galaxy luminosities and the 
gain and loss of galaxies by merging. Here is an example of the former in a spatially-flat cosmological model with 
^mo = 0.25. The expansion time from high redshift is £3 = 2.4 Gyr at redshift z — 3 and to = 15 Gyr now. Consider 
a galaxy observed at z = 3. Suppose the bulk of the stars in this galaxy formed at time tf, and the population then 
aged and faded without significant later star formation. Then if tf <C £3 the ratio of the observed luminosity at z = 3 
to its present luminosity would be (Tinsley, 1972; Worthey, 1994) 

L 3 /L ~ (ioAaf 8 * 4. (75) 

If tf were larger, but still less than £3, this ratio would be larger. If tf were greater than £3 the galaxy would not 
be seen, absent earlier generations of stars. In a more realistic picture significant star formation may be distributed 
over a considerable range of redshifts, and the effect on the typical galaxy luminosity at a given redshift accordingly 
more complicated. Since there are many more galaxies with low luminosities than galaxies with high luminosities, 
one has to know the luminosity evolution quite well for a meaningful comparison of galaxy counts at high and low 
redshifts. The present situation is illustrated by the rather different indications from studies by Phillipps et al. (2000) 
and Totani et al. (2001). 

The understanding of galaxy evolution and the interpretation of galaxy counts will be improved by large samples 
of counts of galaxies as a function of color, apparent magnitude, and redshift. Newman and Davis (2000) point to a 
promising alternative: count galaxies as a function of the internal velocity dispersion that in spirals correlates with 
the dispersion in the dark matter halo. That could eliminate the need to understand the evolution of star populations. 
There is still the issue of evolution of the dark halos by merging and accretion, but that might be reliably modeled 
by numerical simulations within the CDM picture. Either way, with further work galaxy counts may provide an 
important test for dark energy and its evolution (Newman and Davis, 2000; Hutcrer and Turner, 2001; Podariu and 
Ratra, 2001). 



6. The gravitational lensing rate 



The probability of production of multiple images of a quasar or a radio source by gravitational lensing by a 
foreground galaxy, or of strongly lensed images of a galaxy by a foreground cluster of galaxies, adds the relativistic 
expression for the deflection of light to the physics of the homogeneous cosmological model. Fukugita, Futamase, and 
Kasai (1990) and Turner (1990) point out the value of this test: at small Qmo the predicted lensing rate is considerably 
larger in a flat model with A than in an open model with A = (as illustrated in Fig 13.12 in Peebles, 1993). 



Podariu and Ratra (2000) and Waga and Frieman (2000) discuss the redshift-magnitudc relation in the inverse power-law scalar field 
model, and Waga and Frieman (2000) and Ng and Wiltshire (2001) discuss this relation in the massive scalar field model. 

Podariu, Nugent, and Ratra (2001), Weller and Albrecht (2002), Wang and Lovelace (2001), Gerke and Efstathiou (2002), Eriksson 
and Amanullah (2002), and references therein, discuss constraints on cosmological parameters from the proposed SNAP mission. 
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The measurement problem for the analysis of quasar lensing is that quasars that are not lensed are not magnified 
by lensing, making them harder to find and the correction for completeness of detection harder to establish. Present 
estimates (Falco, Kochanek, and Muhoz, 1998; Helbig et ai, 1999) do not seriously constrain f^Mo in an open model, 
and in a flat model (^ko = 0) suggest flyio > 0.36 at 2c This is close to the upper bound in Eq. (p9|) . Earlier 
indications that the lensing rate in a flat model with constant A requires a larger value of Qmo than is suggested 
by galaxy dynamics led Ratra and Quillen (1992) and Waga and Frieman (2000) to investigate the inverse power- 
law potential dark energy scalar field case. They showed this can significantly lower the predicted lensing rate at 
^ko = and small ^mo- The lensing rate still is too uncertain to draw conclusions on this point, but advances in the 
measurement certainly will be followed with interest. 

The main problem in the interpretation of the rate of strong lensing of galaxies by foreground clusters as a cosmo- 
logical test is the sensitivity of the lensing cross section to the mass distribution within the cluster (Wu and Hammer, 
1993); for the present still somewhat uncertain state of the art see Cooray (1999) and references therein. 



7. Dynamics and the mean mass density 

Estimates of the mean mass density from the relation between the mass distribution and the resulting peculiar 
velocities, 81 and from the gravitational deflection of light, probe gravity physics and constrain f^Mo- The former is 
not sensitive to f^KOi ^aoj or the dynamics of the dark energy, the latter only through the angular size distances. 

We begin with the redshift space of observed galaxy angular positions and redshift distances z/Hq in the radial 
direction. The redshift z has a contribution from the radial peculiar velocity, which is a probe of the gravitational 
acceleration produced by the inhomogeneous mass distribution. The two-point correlation function, £ v , in redshift 
space is defined by the probability that a randomly chosen galaxy has a neighbor at distance r\\ along the line of sight 
in redshift space and perpendicular distance r±, 

dP = n(l+£ v (r h r ± ))d ril d 2 r ± , (76) 

where n is the galaxy number density. This is the usual definition of a reduced correlation function. Peculiar velocities 
make the function anisotropic. On small scales the random relative peculiar velocities of the galaxies broaden 
along the line of sight. On large scales the streaming peculiar velocity of convergence to gravitationally growing mass 
concentrations flattens along the line of sight. 82 

At 10 kpc £ hr± ^ 1 Mpc the measured line-of-sight broadening is prominent, and indicates the one-dimensional 
relative velocity dispersion is close to independent of r± at a ~ 300 km s -1 . 83 This is about what would be expected 
if the mass two-and three-point correlation functions were well approximated by the galaxy correlation functions, the 
mass clustering on these scales were close to statistical equilibrium, and the density parameter were in the range of 
Eq. ©. 

We have a check from the motions of the galaxies in and around the Local Group of galaxies, where the absolute 
errors in the measurements of galaxy distances are least. The two largest group members are the Andromeda Nebula 
(M 31) and our Milky Way galaxy. If they contain most of the mass their relative motion is the classical two- 
body problem in Newtonian mechanics (with minor corrections for A, mass accretion at low redshifts, and the tidal 
torques from neighboring galaxies). The two galaxies are separated by 800 kpc and approaching at 110 km s . In 
the minimum mass solution the galaxies have completed just over half an orbit in the cosmological expansion time 
to ~ 10 10 yr. By this argument Kahn and Woltjer (1959) find the sum of masses of the two galaxies has to be an 
order of magnitude larger than what is seen in the luminous parts. An extension to the analysis of the motions and 
distances of the galaxies within 4 Mpc distance from us, and taking account of the gravitational effects of the galaxies 
out to 20 Mpc distance, gives masses quite similar to what Kahn and Woltjer found, and consistent with J^mo i n the 



range of Eq. (59) (Peebles et ai, 2001) 



Early estimates of the mean mass density, by Hubble (1936, p. 189) and Oort (1958), combine the galaxy number density from galaxy 
counts with estimates of galaxy masses from the internal motions of gas and stars. Hubble (1936, p. 180) was quite aware that this 
misses mass between the galaxies, and that the motions of galaxies within clusters suggests there is a lot more intcrgalactic mass (Zwicky, 
1933; Smith, 1936). For a recent review of this subject see Bahcall et. al. (2000). 

This approach grew out of the statistical method introduced by Geller and Peebles (1973); it is derived in its present form in Peebles 
(1980b) and first applied to a serious redshift sample in Davis and Peebles (1983b). These references give the theory for the second 
moment a 2 of in the radial direction — the mean square relative peculiar velocity — in the small-scale stable clustering limit. The 
analysis of the anisotropy of in the linear perturbation theory of large-scale flows (Eq. [p5[) is presented in Kaiser (1987). 

This measurement requires close attention to clusters that contribute little to the mean mass density but a broad and difficult to 
measure tail to the distribution of relative velocities. Details may be traced back through Padilla et al. (2001), Peacock et al. (2001), 
and Landy (2002). 
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We have another check from weak lensing: the shear distortion of images of distant galaxies by the gravitational 
deflection by the inhomogeneous mass distribution. 84 If galaxies trace mass these measurements say the matter density 
parameter measured on scales from about 1 Mpc to 10 Mpc is in the range of Eq. (|59]). It will be interesting to see 
whether these measurements can check the factor of two difference between the relativistic gravitational deflection of 
light and the naive Newtonian deflection angle. 

The redshift space correlation function (Eq. j76|) is measured well enough at hr± ~ 10 Mpc to demonstrate the 
flattening effect, again consistent with £Imo i n the range of Eq. (p39|), if galaxies trace mass. Similar numbers follow 
from galaxies selected as far infrared IRAS sources (Tadros et al, 1999) and from optically selected galaxies (Padilla 
et al., 2001; Peacock et al., 2001). The same physics, applied to estimates of the mean relative peculiar velocity of 
galaxies at separations ~ 10 Mpc, yet again indicates a similar density parameter (Juszkiewicz et al., 2000). 

Other methods of analysis of the distributions of astronomical objects and peculiar velocities smoothed over scales 
£ 10 Mpc give a variety of results for the mass density, some above the range in Eq. ([39]), 85 others towards the 
bottom end of the range (Branchini et al, 2001). The measurement of JImo from large-scale streaming velocities thus 
remains open. But we are impressed by an apparently simple local situation, the peculiar motion of the Local Group 
toward the Virgo cluster of galaxies. This is the nearest known large mass concentration, at distance ~ 20 Mpc. 
Burstein (2000) finds that our virgocentric velocity is v v = 220 km s _1 , indicating SImo — 0.2 (Davis and Peebles, 
1983a, Fig. 1). This leads us to conclude that the weight of the evidence from dynamics on scales ~ 10 Mpc favors 
low r^MOj in the range of Eq. (|59|). 

None of these measurements is precise. But many have been under discussion for a long time and seem to us to 
be reliably understood. Weak lensing is new, but the measurements are checked by several independent groups. The 
result, in our opinion, is a well checked and believable network of evidence that over two decades of well-sampled 
length scales, 100 kpc to 10 Mpc, the apparent value of f^Mo is constant to a factor of three or so, in the range 
0.15 £ f^Mo ~ 0.4. The key point for the purpose of this review is that this result is contrary to what might have been 
expected from biasing, or from a failure of the inverse square law (as will be discussed in test [13]). 



8. The baryon mass fraction in clusters of galaxies 

Abell (1958) made the first useful catalog of the rich clusters considered here and in the next test. A typical value 
of the Abell cluster mass within the Abell radius r a = 1.5/i _1 Mpc is 3 x l0 14 h~ 1 M Q . The cluster masses are reliably 
measured (within Newtonian gravity) from consistent results from the velocities of the galaxies, the pressure of the 
intracluster plasma, and the gravitational deflection of light from background galaxies. 

White (1992) and White et al. (1993) point out that rich clusters likely are large enough to contain a close to fair 
sample of baryons and dark matter, meaning the ratio of baryonic to total mass in a cluster is a good measure of 
f^Bo/^Mo- With f^Bo from the model for light elements (Eqs. p2| ), this gives a measure of the mean mass density. 
The baryon mass fraction in clusters is still under discussion. 8 " We adopt as the most direct and so maybe most 
reliable approach the measurement of the baryonic gas mass fraction of clusters, / gas , through the Sunyaev-Zel'dovich 
microwave decrement caused by Thomson-Compton scattering of cosmic microwave background radiation by the 
intracluster plasma. The Carlstrom et al. (2001) value for / gas gives n M0 ~ 0.25, 87 in the range of Eq. (|59|). This 
test does not directly constrain f^KOi ^A0j or the dynamics of the dark energy. 



9. The cluster mass function 



In the CDM model rich clusters of galaxies grow out of the rare peak upward fluctuations in the primeval Gaussian 
mass distribution. Within this model one can adjust the amplitude of the mass fluctuations to match the abundance 
of clusters at one epoch. In the Einstein-de Sitter model it is difficult to see how this one free adjustment can account 
for the abundance of rich clusters now and at redshifts near unity. 88 



Recent studies include Wilson, Kaiser, and Luppino (2001), Van Waerbeke et al. (2002), Refregier, Rhodes, and Groth (2002), Bacon 
et al. (2002), and Hoekstra, Yee, and Gladders (2002). See Munshi and Wang (2002) and references therein for discussions of how weak 
lensing might probe dark energy. 

The methods and results may be traced through Fisher, Scharf, and Lahav (1994), Sigad et al. (1998), and Branchini et al. (2000). 
See Hradecky et al. (2000), Roussel, Sadat, and Blanchard (2000), Allen, Schmidt, and Fabian (2002), and references therein. 

This assumes Cbo^ 2 = 0.014 from Eqs. j62[). For the full range of values in Eqs. (^j) and (p2]), 0.1 Omo ~ 0.4 at two standard 
deviations. 

Early discussions of this problem include Evrard (1989), Peebles, Daly, and Juszkiewicz (1989), and Oukbir and Blanchard (1992). 
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Most authors now agree that the low density flat ACDM model can give a reasonable fit to the cluster abundances 
as a function of redshift. The constraint on f^Mo from the present cluster abundance still is under discussion, but 
generally is found to be close to Omo ~ 0.3 if galaxies trace mass. 89 The constraint from the evolution of the cluster 
number density also is under discussion. 90 The predicted evolution is slower in a lower density universe, and at given 
f^Mo the evolution is slower in an open model with A = than in a spatially-flat model with A (for the reasons 
discussed in Sec. III.D). Bahcall and Fan (1998) emphasize that we have good evidence for the presence of some 
massive clusters at z <~ 1, and that this is exceedingly difficult to understand in the CDM model in the Einstein-de 
Sitter cosmology (when biasing is adjusted to get a reasonable present number density). Low density models with 
or without A can account for the existence of some massive clusters at high redshift. Distinguishing between the 
predictions of the spatially curved and flat low density cases awaits better measurements. 



10. Biasing and the development of nonlinear mass density fluctuations 

Elements of the physics of cluster formation in test (9) appear in this test of the early stages in the nonlinear growth 
of departures from homogeneity. An initially Gaussian mass distribution becomes skew as low density fluctuations 
start to bottom out and high density fluctuations start to develop into prominent mass peaks. The early signature 
of this nonlinear evolution is the disconnected three-point mass autocorrelation function, (5(x,t)5(y,t)5(z,t)}, where 
S(x,t) = Sp/p is the dimensionless mass contrast. If galaxies are useful mass tracers the galaxy three-point function 
is a good measure of this mass function. 

The form for the mass three-point function, for Gaussian initial conditions at high redshift, in lowest nonzero order 
in perturbation theory, is worked out in Fry (1984), and Fry (1994) makes the point that measurements of the galaxy 
three-point function test how well galaxies trace mass. 91 There are now two sets of measurements of the galaxy 
three-point function on scales ~ 10 to 20 Mpc, where the density fluctuations are not far from Gaussian. One uses 
infrared-selected IRAS galaxies, 92 the other optically-selected galaxies (Verde et al, 2002). The latter is consistent 
with the perturbative computation of the mass three-point function for Gaussian initial conditions. The former says 
infrared-selected galaxies are adequate mass tracers apart from the densest regions, which IRAS galaxies avoid. That 
has a simple interpretation in astrophysics: galaxies in dense regions tend to be swept clear of the gas and dust that 
make galaxies luminous in the infrared. 

This test gives evidence of consistency of three ideas: galaxies are useful mass tracers on scales ~ 10 Mpc, the initial 
conditions are close to Gaussian, and conventional gravity physics gives an adequate description of this aspect of the 
growth of structure. It is in principle sensitive to I^ao, through the suppression of the growth of small departures 
from homogeneity at low redshift, but the effect is small. 



11. The anisotropy of the cosmic microwave background radiation 

The wonderfully successful CDM prediction of the power spectrum of the angular distribution of the temperature 
of the 3 K cosmic microwave background radiation has converted many of the remaining skeptics in the cosmology 
community to the belief that the CDM model likely does capture important elements of reality. 

Efstathiou (2002) provides a useful measure of the information in the present measurements 93 : the fit to the 
CDM model significantly constrains three linear combinations of the free parameters. We shall present three sets 
of considerations that roughly follow Efstathiou's constraints. We begin with reviews of the standard measure of 
the temperature anisotropy and of the conditions at redshift z ~ 1000 that are thought to produce the observed 
anisotropy. 

The 3 K cosmic microwave background temperature T{9, <p) as a function of position in the sky usually is expressed 



For recent discussions see Pierpaoli, Scott, and White (2001), Seljak (2001), Viana, Nichol, and Liddle (2002), Ikebe et al. (2002), 
Bahcall et al. (2002), and references therein. Wang and Stcinhardt (1998) consider this test in the context of the XCDM paramctrization; 
to our knowledge it has not been studied in the scalar field dark energy case. 

Examples include Blanchard et al. (2000), and Borgani et al. (2001). 

Other notable contributions to the development of this point include Bernardeau and SchaefTer (1992), Fry and Gaztanaga (1993), and 
Hivon et al. (1995). 

Two sub-samples of IRAS galaxies are analyzed by Scoccimarro et al. (2001) and by Fcldman et al. (2001). 

Recent measurements are presented in Lee et al. (2001), Ncttcrfield et al. (2002), Halverson et al. (2002), Miller et al. (2002a), Coble 
et al. (2001), Scott et al. (2002), and Mason et al. (2002). 
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as an expansion in spherical harmonics, 

ST(0, 4>) = T(0, 0) - (T) = £ opy,"^, 0). (77) 

? . m 

The square of 8T averaged over the sky is 

< (5r2 ) = ^E( 2Z + 1 )(i a ™i) 2 ' ( 78 ) 
i 

where |a™| 2 is statistically independent of m. This may be rewritten as 

(5T 2 )=Y / 1 1 STt, 5T? = l V±p±{\aT?). (79) 
i 

Since X ^ s close to J d\nl, 8T? is the variance of the temperature per logarithmic interval of I. A measure of the 
angular scale belonging to the multipole index I is that the minimum distance between zeros of the spherical harmonic 
Y™, in longitude or latitude, is 6 — n/l, except close to the poles, where Yj m approaches zero. 94 

Now let us consider the main elements of the physics that determines the 3 K cosmic microwave background 
anisotropy. 95 At redshift Zd cc ~ 1000 the temperature reaches the critical value at which the primeval plasma 
combines to atomic hydrogen (and slightly earlier to neutral helium). This removes the coupling between baryons and 
radiation by Thomson scattering, leaving the radiation to propagate nearly freely (apart from residual gravitational 
perturbations) . Ratios of mass densities near the epoch z^cc when matter and radiation decouple are worth noting. 
At redshift z cq — 2.4 x 10 4 f2Mo^ 2 the mass density in matter — including the baryonic and nonbaryonic components 
— is equal to the relativistic mass density in radiation and neutrinos assumed to have low masses. At decoupling the 
ratio of mass densities is 

^dec2 = £eq_ _ _ ^ (8Q) 

PR(Zdcc) ^dcc 

at the central values of the parameters in Eqs. (|J) and (|59"1). The ratio of mass densities in baryons and in thermal 
cosmic microwave background radiation — not counting neutrinos — is 

Pb{z6cc) _ 4 x 10 4 ^bo^_ _ 5 , gl s 



PCBR^dec) 

That is, the baryons and radiation decouple just as the expansion rate has become dominated by nonrelativistic 
matter and the baryons are starting to lower the velocity of sound in the coupled baryon-radiation fluid (presenting 
us with still more cosmic coincidences). 

The acoustic peaks in the spectrum of angular fluctuations of the 3 K cosmic microwave background radiation come 
from the Fourier modes of the coupled baryon-radiation fluid that have reached maximum or minimum amplitude 
at decoupling. Since all Fourier components start at zero amplitude at high redshift — in the growing density 
perturbation mode — this condition is 

/ kc s dt/a ~ nn/2, (82) 
Jo 

where c s is the velocity of sound in the baryon-radiation fluid. Before decoupling the mass density in radiation is 
greater than that of the baryons, so the velocity of sound is close to c/ V3. The proper wavelength at the first acoustic 
peak thus is 

Apeak ~ ^dec h 1 Qj v iq ■ (83) 



A more careful analysis distinguishes averages across the sky from ensemble averages. By historical accident the conventional normal- 
ization replaces 21 + 1 with 2(1 + 1) in Eq. (|79[). Kosowsky (2002) reviews the physics of the polarization of the radiation. 

The physics is worked out in Peebles anaYu (1970) and Peebles (1982). Important analytic considerations are in Sunyaev and 
Zel'dovich (1970). The relation of the cosmic microwave background anisotropy to the cosmological parameters is explored in many 
papers; examples of the development of ideas include Bond (1988), Bond et al. (1994), Hu and Sugiyama (1996), Ratra et al. (1997, 
1999), Zaldarriaga, Spergel, and Scljak (1997), and references therein. 
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The parameter dependence comes from Eq. (^36|). The observed angle subtended by A pca k is set by the angular size 
distance r computed from z cq to the present (Eq. fn]] ). If f2xo = or fl A0 = the angular size distance is 

H a r ~ 2ft-™, m = 1 if Q A0 = 0, m ~ 0.4 if ft K0 = 0. (84) 

If A = this expression is analytic at large z cq . The expression for Qko = is a reasonable approximation to the 
numerical solution. So the angular scale of the peak varies with the matter density parameter as 

1/2 

fpeak ^ ^dec Apeak/ao^ oc ft^ if ft A o = 0, (85) 
cx n^- 1 if n K0 = 0. (86) 

The key point from these considerations is that the angle defined by the first peak in the fluctuation power spectrum 
is sensitive to J7mo if A = (Eq. [pS]] ), but not if £1ko = (Eq. |8(|). 96 We have ignored the sensitivity of Zdcc and idee 
to Qmo, but the effect is weak. More detailed computations, which are needed for a precise comparison with the data, 

— 1/2 

show that the CDM model predicts that the first and largest peak of STi appears at multipole index Z pca k ~ 220£1 MO 
if A = 0, and at Z pcak ~ 220 if Q K0 = and 0.1 < ft M0 ~ I- 97 

The measured spectrum 98 peaks at STi ~ 80 /zK at / ~ 200, thus requiring small space curvature in the CDM 
model. This is the first of Efstathiou's constraints. Because of the geometric degeneracy this measurement does not 
yet seriously constrain Omo if ^ko = 0. 

The second constraint comes from the spectrum of temperature fluctuations on large scales, Z 30, where pressure 
gradient forces never were very important. Under the scale-invariant initial conditions discussed in Sec. III.C the 
Einstein-de Sitter model predicts 6Ti is nearly independent of / on large scales. A spatially-flat model with J7mo ~ 0.3, 
predicts <5T; decreases slowly with increasing I at small 1. 99 The measured spectrum is close to flat at 6Ti ~ 30 fiK, 
but not well enough constrained for a useful measure of the parameters Qmq and SIao. 100 Because of the simplicity of 
the physics on large angular scales, this provides the most direct and so perhaps most reliable normalization of the 
CDM model power spectrum (that is, the parameter A in Eqs. jl0| and p|). 

The third constraint is the baryon mass density. It affects the speed of sound c s (Eq. ||32fl) in the baryon-radiation 
fluid prior to decoupling, and the mean free path for the radiation at z ~ Zdec- These in turn affect the predicted 
sequence of acoustic peaks (see, e.g., Hu and Sugiyama, 1996). The detected peaks are consistent with a value 
for the baryon density parameter J7bo in a range that includes what is derived from the light elements abundances 
(Eqs. f32"| ). 101 This impressive check may be much improved by the measurements of STi in progress. 

The measurements of STi are consistent with a near scale- invariant power spectrum (Eq. |4l|] with n ~ 1) with 
negligible contribution from gravity wave or isocurvature fluctuations (Sec. III.C.l). The 3 K cosmic microwave 
background temperature fluctuations show no departure from a Gaussian random process. 102 This agrees with the 



This "geometrical degeneracy" is discussed by Efstathiou and Bond (1999). Marriage (2002) presents a closer analysis of the effect. 
Sugiyama and Gouda (1992), Kamionkowski, Spergel, and Sugiyama (1994b), and Kamionkowski et al. (1994a) are early discussions of 
the cosmic microwave background anisotropy in an open model. 

Brax et al. (2000) and Baccigalupi et al. (2000) compute the angular spectrum of the cosmic microwave background anisotropy in the 
dark energy scalar field model. Doran et al. (2001) discuss the angular scale of the peaks in this case, and Corasaniti and Copeland 
(2002), Baccigalupi et al. (2002), and references therein, compare model predictions and observations — it is too early to draw profound 
conclusions about model viability, and new data are eagerly anticipated. Wasserman (2002) notes that the cosmic microwave background 
anisotropy data could help discriminate between different dark energy scalar field models whose predictions do not differ significantly 
at low redshift. 

For analyses see Knox and Page (2000), Podariu et al. (2001), Wang, Tegmark, and Zaldarriaga (2002), Durrer, Novosyadlyj, and 
Apunevych (2002), Miller et al. (2002b), and references therein. 

The physics was first demonstrated by Sachs and Wolfe (1967) and applied in the modern context by Peebles (1982). The intermediate 
Sachs- Wolfe effect that applies if the universe is not Einstein-de Sitter is shown in Eq. (93.26) in Peebles (1980). This part of the 
Sachs- Wolfe effect receives a contribution from the low redshift matter distribution, so cross-correlating the observed large-scale cosmic 
microwave background anisotropy with the low redshift matter distribution could provide another test of the world model (Boughn and 
Crittenden, 2001, and references therein). 

See, e.g., Gorski et al. (1998). This ignores the "low" value of the cosmological quadrupolc (I = 2) moment, whose value depends on 
the model used to remove foreground Galactic emission (see, e.g., Kogut et al., 1996). Contamination due to non-cosmic microwave 
background emission is an issue for some of the anisotropy data sets (see, e.g., de Oliveira-Costa et al., 1998; Hamilton and Ganga, 2001; 
Mukherjee et al., 2002, and references therein). Other issues that need care in such analyses include accounting for the uncertainty in 
the calibration of the experiment (see, e.g., Ganga et al, 1997; Bridle et al., 2002), and accounting for the shape of the antenna pattern 
(see, e.g., Wu et al., 2001a; Souradeep and Ratra, 2001; Fosalba, Dore, and Bouchct, 2002). 

The Qbo^ 2 values estimated from the cosmic microwave background anisotropy measured by Netterfield et al. (2002), P rvke et al. 
(2002), and Stompor et al. (2001), are more consistent with the higher, deuterium based, Buries et al. (2001) range in Eqs. (k32|). 

Colley, Gott, and Park (1996) present an early discussion of the situation on large angular scales; more recent discussions are in 
Mukherjee, Hobson, and Lasenby (2000), Phillips and Kogut (2001), and Komatsu et al. (2002). Degree and sub-degree angular scale 
anisotropy data are studied in Park et al. (2001), Wu et al. (2001b), Shandarin et al. (2002), and Polenta et al. (2002). 



40 



picture in test (10) for the nonlinear growth of structure out of Gaussian initial mass density fluctuations. 

The interpretation of the cosmic microwave background temperature anisotropy measurements assumes and tests 
general relativity and the CDM model. One can write down other models for structure formation that put the peak 
of STi at about the observed angular scale — an example is Hu and Peebles (2000) — but we have seen none so far 
that seem likely to fit the present measurements of 5Ti. Delayed recombination of the primeval plasma in an low 
density A — CDM model can shift the peak of STi to the observed scale. 103 The physics is valid, but the scenario 
is speculative and arguably quite improbable. On the other hand, we cannot be sure a fix of the challenges to CDM 
reviewed in Sec. IV. A. 2 will not affect our assessments of such issues, and hence of this cosmological test. 



12. The mass autocorrelation function and nonbaryonic matter 

If the bulk of the nonrelativistic matter, with density parameter J7mo ~ 0.25, were baryonic, then under adiabatic 
initial conditions the most immediate problem would be the strong dissipation of primeval mass density fluctuations 
on the scale of galaxies by diffusion of radiation through the baryons at redshifts near decoupling. 104 Galaxies could 
form by fragmentation of the first generation of protocluster "pancakes," as Zel'dovich (1978) proposed, but this 
picture is seriously challenged by the evidence that the galaxies formed before clusters of galaxies. 105 In a baryonic 
dark matter model we could accommodate the observations of galaxies already present at z ~ 3 by tilting the primeval 
mass fluctuation spectrum to favor large fluctuations on small scales, but that would mess up the cosmic microwave 
background anisotropy. The search for isocurvature initial conditions that might fit both in a baryonic dark matter 
model has borne no fruit so far (Peebles, 1987). 

The most important point of this test is the great difficulty of reconciling the power spectra of matter and radiation 
without the postulate of nonbaryonic dark matter. The CDM model allows hierarchical growth of structure, from 
galaxies up, which is what seems to be observed, because the nonbaryonic dark matter interacts with baryons and 
radiation only by gravity; the dark matter distribution is not smoothed by the dissipation of density fluctuations in 
the baryon-radiation fluid at redshifts z ^ z cq . 

As discussed in Sec. III.D, in the CDM model the small scale part of the dark matter power spectrum bends from 
the primeval scale-invariant form P(k) oc k to P(k) oc k~ 3 , and the characteristic length at the break scales inversely 
with r^Mu (Eq. p^f ). Evidence of such a break in the galaxy power spectrum P g (k) has been known for more than a 
decade 106 ; it is consistent with a value of f^Mo m the range of Eq. (p9|). 



13. The gravitational inverse square law 

The inverse square law for gravity determines the relation between the mass distribution and the gravitationally- 
driven peculiar velocities that enter estimates of the matter density parameter f^Mo- The peculiar velocities also 
figure in the evolution of the mass distribution, and hence the relation between the present mass fluctuation spectrum 
and the spectrum of cosmic microwave background temperature fluctuations imprinted at redshift z ~ 1000. We are 
starting to see demanding tests of both aspects of the inverse square law. 

We have a reasonably well checked set of measurements of the apparent value of Omo on scales ranging from 100 kpc 
to 10 Mpc (as reviewed under test [7]). Most agree with a constant value of the apparent ^mo, within a factor of 
three or so. This is not the precision one would like, but the subject has been under discussion for a long time, and, 
we believe, is now pretty reliably understood, within the factor of three or so. If galaxies were biased tracers of mass 
one might have expected to have seen that Omo increases with increasing length scale, as the increasing scale includes 
the outer parts of extended massive halos. Maybe that is masked by a gravitational force law that decreases more 
rapidly than the inverse square law at large distance. But the much more straightforward reading is that the slow 



The model in Peebles, Seager, and Hu (2000) assumes stellar ionizing radiation at z ~ 1000 produces recombination Lyman a photons. 
These resonance photons promote photoionization from the n = 2 level of atomic hydrogen. That allows delayed recombination with a 
rapid transition to neutral atomic hydrogen, as required to get the shape of ST t about right. 
Early analyses of this effect are in Peebles (1965), and Silk (1967, 1968). 

For example, our Milky Way galaxy is in the Local Group, which seems to be just forming, because the time for a group member to 
cross the Local Group is comparable to the Hubble time. The Local Group is on the outskirts of the concentration of galaxies around 
the Virgo cluster. We and neighboring galaxies are moving away from the cluster, but at about 80 percent of the mean Hubble flow, 
as if the local mass concentration were slowing the local expansion. That is, our galaxy, which is old, is starting to cluster with other 
galaxies, in a "bottom up" hierarchical growth of structure, as opposed to the "top down" evolution of the pancake picture. 

The first good evidence is discussed in Efstathiou et al. (1990); for recent examples see Sutherland et al. (1999), Percival et al. (2001), 
and Dodelson et al. (2002). 
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variation of f^Mo sampled over two orders of magnitude in length scale agrees with the evidence from tests (7) to (10) 
that galaxies are useful mass tracers, and that the inverse square law therefore is a useful approximation on these 
scales. 

The toy model in Eq. (|57]) illustrates how a failure of the inverse square law would affect the evolution of the shape 
of the mass fluctuation power spectrum P(k,t) as a function of the comoving wavenumber k, in linear perturbation 
theory. This is tested by the measurements of the mass and cosmic microwave background temperature fluctuation 
power spectra. The galaxy power spectrum P g (k) varies with wavenumber at k ~ O.lh Mpc -1 about as expected 
under the assumptions that the mass distribution grew by gravity out of adiabatic scale-invariant initial conditions, 
the mass is dominated by dark matter that does not suffer radiation drag at high redshift, the galaxies are useful 
tracers of the present mass distribution, the matter density parameter is JImo ~ 0.3, and, of course, the evolution is 
adequately described by conventional physics (Hamilton and Tegmark, 2002, and references therein). If the inverse 
square law were significantly wrong at k ~ O.lh Mpc , the near scale-invariant form would have to be an accidental 
effect of some failure in this rather long list of assumptions. This seems unlikely, but a check certainly is desirable. 
We have one, from the cosmic microwave background anisotropy measurements. They also are consistent with near 
scale-invariant initial conditions applied at redshift z ~ 1000. This preliminary check on the effect of the gravitational 
inverse square law applied on cosmological length scales and back to redshift z ~ 1000 will be improved by better 
understanding of the effect on STi of primeval tensor perturbations to spacetime, and of the dynamical response of 
the dark energy distribution to the large-scale mass distribution. 

Another aspect of this check is the comparison of values of the large-scale rms fluctuations in the present distributions 
of mass and the cosmic microwave background radiation. The latter is largely set at decoupling, after which the former 
grows by a factor of about 10 3 to the present epoch, in the standard relativistic cosmological model. If space curvature 
is negligible the growth factor agrees with the observations to about 30%, assuming galaxies trace mass. In a low 
density universe with A = the standard model requires that mass is more smoothly distributed than galaxies, 
5N/N r~j 3SM/M, or that the gravitational growth factor since decoupling is a factor of three off the predicted 
factor ~ 1000; this factor of three is about as large a deviation from unity as is viable. We are not proposing this 
interpretation of the data, rather we are impressed by the modest size of the allowed adjustment to the inverse square 
law. 



C. The state of the cosmological tests 

Precision cosmology is not very interesting if it is based on faulty physics or astronomy. That is why we have 
emphasized the tests of the standard gravity physics and structure formation model, and the checks of consistency of 
measures based on different aspects of the astronomy. 

There are now five main lines of evidence that significantly constrain the value of ^mo to the range of Eq. (|5{]): the 
redshift-magnitude relation (test [4]), gravitational dynamics and weak lensing (test [7]), the baryon mass fraction 
in clusters of galaxies (test [8]), the abundance of clusters as a function of mass and redshift (test [9]), and the 
large-scale galaxy distribution (test [12]). There are indications for larger values of f^Mo, from analyses of the rate 
of strong lensing of quasars by foreground galaxies (test [6]) and some analyses of large-scale flows (test [7]), though 
we know of no well-developed line of evidence that points to the Einstein-de Sitter value Qmo — 1. Each of these 
measures of f^Mo may suffer from systematic errors: we must bear in mind the tantalus principle mentioned in Sec. 
LA, and we have to remember that the interpretations could be corrupted by a failure of standard physics. But the 
general pattern of results from a considerable variety of independent approaches seems so close to consistent as to 
be persuasive. Thus we conclude that there is a well-checked scientific case for the proposition that the measures of 
the mean mass density of matter in forms capable of clustering are physically meaningful, and that the mass density 
parameter almost certainly is in the range 0.15 £ £Imq ~ 0.4. 

In the standard cosmology the masses of the galaxies are dominated by dark matter, with mass density parameter 
Admo ~ 0.2, that is not baryonic (or acts that way). We do not have the direct evidence of a laboratory detection; 
this is based on two indirect lines of argument. First, the successful model for the origin of the light elements (test [2]) 
requires baryon density fieo ~ 0.05. It is difficult to see how to reconcile a mass density this small with the mass 
estimates from dynamics and lensing; the hypothesis that f^Mo is dominated by matter that is not baryonic allows 
us to account for the difference. Second, the nonbaryonic matter allows us to reconcile the theory of the anisotropy 
of the cosmic microwave background radiation with the distributions of galaxies and groups and clusters of galaxies, 
and the presence of galaxies at z ~ 3 (tests [11] and [12]). This interpretation requires a value for fieo that is in line 
with test (2). The consistency is impressive. But the case is not yet as convincing as the larger network of evidence 
that £Imo is well below unity. 

The subject of this review is Einstein's cosmological constant A, or its equivalent in dark energy. The evidence for 
detection of A by the redshift-magnitude relation for type la supernovae is checked by the angular distribution of the 
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3 K cosmic microwave background temperature together with the constraints on Omo- This certainly makes a serious 
case for dark energy. But we keep accounts by the number of significant independent checks, and by this reckoning 
the case is not yet as strong as for nonbaryonic dark matter. 

V. CONCLUDING REMARKS 

We cannot demonstrate that there is not some other physics, applied to some other cosmology, that equally well 
agrees with the cosmological tests. The same applies to the whole enterprise of physical science, of course. Parts of 
physics are so densely checked as to be quite convincing approximations to the way the world really is. The web of 
tests is much less dense in cosmology, but, we have tried to demonstrate, by no means negligible, and growing tighter. 

A decade ago there was not much discussion of how to test general relativity theory on the scales of cosmology. 
That was in part because the theory seems so logically compelling, and certainly in part also because there was not 
much evidence to work with. The empirical situation is much better now. We mentioned two tests, of the relativistic 
active gravitational mass density and of the gravitational inverse square law. The consistency of constraints on f2 M0 
from dynamics and the redshift-magnitude relation adds a test of the effect of space curvature on the expansion rate. 
These tests are developing; they will be greatly improved by work in progress. There certainly may be surprises in 
the gravity physics of cosmology at redshifts z £ 10 10 , but it already is clear that if so the surprises are subtly hidden. 

A decade ago it was not at all clear which direction the theory of large-scale structure formation would take. Now 
the simple CDM model has proved to be successful enough that there is good reason to expect the standard model ten 
years from now will resemble CDM. We have listed challenges to this structure formation model. Some may well be 
only apparent, a result of the complications of interpreting the theory and observations. Others may prove to be real 
and drive adjustments of the model. Fixes certainly will include one element from the ideas of structure formation 
that were under debate a decade ago: explosions that rearrange matter in ways that are difficult to compute. Fixes 
may also include primeval isocurvature departures from homogeneity, as in spacetime curvature fluctuations frozen 
in during inflation, and maybe in new cosmic fields. It would not be surprising if cosmic field defects, that have such 
a good pedigree from particle physics, also find a role in structure formation. And a central point to bear in mind is 
that fixes, which do not seem unlikely, could affect the cosmological tests. 

A decade ago we had significant results from the cosmological tests. For example, estimates of the product H to 
suggested we might need positive A, though the precision was not quite adequate for a convincing case. That still is 
so; the community will be watching for further advances. We had pretty good constraints on f^Bo from the theory of 
the origin of the light elements. The abundance measurements are improving; an important recent development is the 
detection of deuterium in gas clouds at redshifts z <~ 3. Ten years ago we had useful estimates of masses from peculiar 
motions on relatively small scales, but more mixed messages from larger scales. The story seems more coherent now, 
but there still is room for improved consistency. We had a case for nonbaryonic dark matter, from the constraint on 
fieo and from the CDM model for structure formation. The case is tighter now, most notably due to the successful 
fit of the CDM model prediction to the measurements of the power spectrum of the anisotropy of the 3 K cosmic 
microwave background radiation. In 1990 there were believable observations of galaxies identified as radio sources at 
z ~ 3. Now the distributions of galaxies and the intergalactic medium at z ~ 3 are mapped out in impressive detail, 
and we are seeing the development of a semi-empirical picture of galaxy formation and evolution. Perhaps that will 
lead us back to the old dream of using galaxies as markers for the cosmological tests. 

Until recently it made sense to consider the constraints on one or two of the parameters of the cosmology while 
holding all the rest at "reasonable" values. That helps us understand what the measurements are probing; it is 
the path we have followed in Sec. IV. B. But the modern and very sensible trend is to consider joint fits of large 
numbers of parameters to the full suite of observations. 107 This includes a measure of the biasing or antibiasing of 
the distribution of galaxies relative to mass, rather than our qualitative argument that one usefully approximates 
the other. In a fully satisfactory cosmological test the parameter set will also include parametrized departures from 
standard physics extrapolated to the scales of cosmology. 

Community responses to advances in empirical evidence are not always close to linear. The popularity of the 
Einstcin-dc Sitter model continued longer into the 1990s than seems logical to us, and the switch to the now standard 
ACDM cosmological model — with flat space sections, nonbaryonic cold dark matter, and dark energy — arguably 
is more abrupt than warranted by the advances in the evidence. Our review leads us to conclude that there is now 
a good scientific case that the matter density parameter is J7mo ~ 0.25, and a pretty good case that about three 
quarters of that is not baryonic. The cases for dark energy and for the ACDM model are significant, too, though 



Examples are Cole et al. (1997), Jenkins et al. (1998), Lineweavcr (2001), and Percival et al. (2002). 
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beclouded by observational issues of whether we have an adequate picture of structure formation. But we expect the 
rapid advances in the observations of structure formation will soon dissipate these clouds, and, considering the record, 
likely reveal new clouds on the standard model for cosmology a decade from now. 

A decade ago the high energy physics community had a well-defined challenge: show why the dark energy density 
vanishes. Now there seems to be a new challenge and clue: show why the dark energy density is exceedingly small but 
not zero. The present state of ideas can be compared to the state of research on structure formation a decade ago: in 
both situations there are many lines of thought but not a clear picture of which is the best direction to take. The big 
difference is that a decade ago we could be reasonably sure that observations in progress would guide us to a better 
understanding of how structure formed. Untangling the physics of dark matter and dark energy and their role in 
gravity physics is a much more subtle challenge, but, we may hope, will also be guided by advances in the exploration 
of the phenomenology. Perhaps in another ten years that will include detection of evolution of the dark energy, and 
maybe detection of the gravitational response of the dark energy distribution to the large-scale mass distribution. 
There may be three unrelated phenomena to deal with: dark energy, dark matter, and a vanishing sum of zero-point 
energies and whatever goes with them. Or the phenomena may be related. Because our only evidence of dark matter 
and dark energy is from their gravity it is a natural and efficient first step to suppose their properties are as simple as 
allowed by the phenomenology But it makes sense to watch for hints of more complex physics within the dark sector. 

The past eight decades have seen steady advances in the technology of application of the cosmological tests, from 
telescopes to computers; advances in the theoretical concepts underlying the tests; and progress through the learning 
curves on how to apply the concepts and technology. We see the results: the basis for cosmology is much firmer than 
a decade ago. And the basis surely will be a lot more solid a decade from now. 

Einstein's cosmological constant, and the modern variant, dark energy, have figured in a broad range of topics in 
physics and astronomy that have been under discussion, in at least some circles, much of the time for the past eight 
decades. Many of these issues undoubtedly have been discovered more than once. But in our experience such ideas 
tend to persist for a long time at low visibility and sometimes low fidelity Thus the community has been very well 
prepared for the present evidence for detection of dark energy. And for the same reason we believe that dark energy, 
whether constant, or rolling toward zero, or maybe even increasing, still will be an active topic of research, in at 
least some circles, a decade from now, whatever the outcome of the present work on the cosmological tests. Though 
this much is clear, we see no basis for a prediction of whether the standard cosmology a decade from now will be a 
straightforward elaboration of ACDM, or whether there will be more substantial changes of direction. 
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VI. APPENDIX: RECENT DARK ENERGY SCALAR FIELD RESEARCH 

Many dark energy models are characterized by attractor or tracker behavior. The goal also is to design the model 
so that the field energy density is subdominant at high rcdshift, where it is not wanted, and dominant at low redshift, 
where that is what seems to be observed. 

In the simplest such scalar field models the action has a conventionally normalized scalar field kinetic and spatial 
gradient term, and the real scalar field is coupled only to itself and gravity. Then the scalar field part of the model 
is fully characterized by the scalar field potential (along with some broad constraints on the initial conditions for the 
field, if the attractor behavior is realized). The inverse power law potential model is described in Sees. II. C and III.E, 
and used in some of the cosmological tests in Sec. IV. Here we list other scalar field potentials now under discussion 
for these minimal dark energy models, modifications of the kinetic part of the action, possible guidance from high 
energy particle physics ideas, and constraints on these ideas from cosmological observations in the context of dark 
energy models. 

To those of us not active in this field the models may seem baroque in their complexity, but that may be the way 
Nature is. And as we accumulate more and better data it will be possible to test and constrain these models. 

This is an active area of research, with frequent introduction of new models, so our discussion must be somewhat 
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sketchy. 108 We limit discussion to a cosmological model with space sections that are flat thanks to the presence of 
the dark energy density. 

As mentioned in Sec III.E, the simplest exponential potential scalar field model is unacceptable because it cannot 
produce the wanted transition from sub-dominant to dominant energy density. Ratra and Peebles (1988) consider 
more complex potentials, such as powers of linear combinations of exponentials of the field Related models are 
under active investigation. These include potentials that are powers of cosh($) and/or sinh($). 109 Sahni and Wang 
(2000) present a detailed discussion of a specific example, V($) cx (coshA$ — l) p , where A and p are constants, 
emphasizing that this potential interpolates from V oc e~ pA * to V cx (A$) 2p as |A$| increases, thus preserving some 
of the desirable properties of the simplest exponential potential case. De la Macorra and Piccinelli (2000) consider 
potentials that are exponentials of more complicated functi ons of $, such as $ 2 and e*. Skordis and Albrecht 
(2002) discuss a model with V($) cx [1 + ($ - A) 2 ]exp[-$y/q/2], while Dodelson, Kaplinghat, and Stewart (2000) 
study a potential V r ($) cx [1 + Asin(_B$)]exp[— $-y/g/2], and Ng, Nunes, and Rosati (2001) consider a model with 
V($) cx § A exp{B$ c ], a simple example of a class of supergravity-inspired potentials studied by Copeland, Nunes, 
and Rosati (2000). Here A, B, and C are parameters. Steinhardt et al. (1999) consider more complicated functions 
of inverse powers of $, such as V r ($) oc cxp(l/$) and linear combinations of inverse powers of $. Brax and Martin 
(2000) consider a supergravity-inspired generalization with V(<&) = K$ _Q e* / 2 . 110 

It will be quite a challenge to select from this wide range of functional forms for the potential those that have 
particular theoretical merit and some chance of being observationally acceptable. 

Following Dolgov (1983), there have been discussions of non-minimally coupled dark energy scalar field models. 111 
These Jordan-Brans-Dicke type models have an explicit coupling between the Ricci scalar and a function of the scalar 
field. This causes the effective gravitational "constant" G (in units where masses are constant) to vary with time, 
which may be an observationally interesting effect and a useful constraint. 112 

In yet another approach, people have considered modifying the form of the dark energy scalar field kinetic and 
spatial gradient term in the action. 113 

We discussed in Sec. III.E the idea that at the end of inflation the dark energy scalar field potential might patch 
on to the part of the scalar field potential responsible for inflation. In an inverse power law model for the dark 
energy potential function this requires an abrupt drop in V($>) at the end of inflation. More sophisticated models, 
now dubbed quintessential inflation, attempt to smooth out this drop by constructing scalar field potentials that 
interpolate smoothly between the part responsible for inflation and the low redshift, dark energy, part. These models 
assume either minimally or non-minimally (Jordan-Brans-Dicke) coupled scalar fields. 114 

The dark energy scalar field models we have reviewed here are meant to be classical, effective, descriptions of what 
might come out of a more fundamental quantum mechanical theory. The effective dark energy scalar field is coupled 
to itself and gravity, and is supposed to be coupled to the other fields in the universe only by gravity. This might 
be what Nature chooses, but we lack an understanding of why the coupling of dark and ordinary fields that are 
allowed by the symmetries are not present, or have coupling strengths that are well below what might be expected 
by naive dimensional analysis (e.g., Kolda and Lyth, 1999). A satisfactory solution remains elusive. Perhaps this is 
not unexpected, because it likely requires a proper understanding of how to reconcile general relativity and quantum 
mechanics. 

We turn now to scalar field dark energy models that arguably are inspired by particle physics. Inverse power 
law scalar field potentials are generated non-perturbatively in models of dynamical supersymmetry breaking. In 



08 In particular, we do not discuss variable mass term models, complex or multiple scalar field models, or repulsive matter. We also omit 
non-scalar field aspects of brane models, Kaluza-Klein models, bimctric theories of gravitation, quantum mechanical running of the 
cosmological constant, the Chaplygin gas, and the supcrstring tachyon. 

09 For examples sec Chimento and Jakubi (1996), Starobinsky (1998), Kruger and Norbury (2000), Di Pietro and Demaret (2001), 
Urena-Lopez and Matos (2000), Gonzalez-Diaz (2000), and Johri (2001). 

10 For still more examples see Green and Lidscy (2000), Barreiro, Copeland, and Nunes (2000), Rubano and Scudcllaro (2002), Sen and 
Sethi (2002), and references therein. 

11 See Uzan (1999), Chiba (1999), Amcndola (1999, 2000), Pcrrotta, Baccigalupi, and Matarrese (2000), Bartolo and Pictroni (2000), 
Bertolami and Martins (2000), Fujii (2000), Faraoni (2000), Baccigalupi et al. (2000), Chen, Scherrer, and Steigman (2001), and 
references therein. 

12 Nucleosynthesis constraints on these and related models are discussed by Arai, Hashimoto, and Fukui (1987), Etoh et al. (1997), 
Perrotta et al. (2000), Chen et al. (2001), and Yahiro et al. (2002). 

13 See Fujii and Nishioka (1990), Chiba, Okabe, and Yamaguchi (2000), Armendariz-Picon, Mukhanov, and Steinhardt (2001), Hebcckcr 
and Wetterich (2001), and references therein. 

14 Early work includes Frewin and Lidsey (1993), Spokoiny (1993), Joyce and Prokopec (1998), and Peebles and Vilenkin (1999); more 
recent discussions are given by Kaganovich (2001), Huey and Lidsey (2001), Majumdar (2001), Sahni, Sami, and Souradccp (2002), and 
Dimopoulos and Valle (2001). 
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supersymmetric non-Abelian gauge theories, the resulting scalar field potential may be viewed as being generated 
by instantons, the potential being proportional to a power of e^ 1 ^ 9 where g(§) is the gauge coupling constant 
which evolves logarithmically with the scalar field through the renormalization group equation. Depending on the 
parameters of the model, an inverse power-law scalar field potential can result. This mechanism may be embedded 
in supergravity and superstring/M theory models. 115 This has not yet led to a model that might be compared to the 
observations. 116 

In the model considered by Weiss (1987), Frieman et al. (1995), and others the dark energy field potential is of 
the form V(&) = M 4 [cos($//) + 1], where M and / are mass scales and the mass of the inhomogeneous scalar field 
fluctuation ~ A/ 2 // is on the order of the present value of the Hubble parameter. For discussions of how this model 
might be more firmly placed on a particle physics foundation see Kim (2000), Choi (1999), Nomura, Watari, and 
Yanagida (2000), and Barr and Seckel (2001). 

There has been much recent interest in the idea of inflation in the brane scenario. Dvali and Tye (1999) note that 
the potential of the scalar field which describes the relative separation between branes can be of a form that leads 
to inflation, and will include some inverse power-law scalar field terms. 117 It will be interesting to learn whether 
these considerations can lead to a viable dark energy scalar field model. Brane models allow for a number of other 
possibilities for dark energy scalar fields, 118 but it is too early to decide whether any of these options give rise to 
observationally acceptable dark energy scalar field models. 

Building on earlier work 119 , Hellerman, Kaloper, and Susskind (2001), and Fischler et al. (2001) note that dark 
energy scalar field cosmological models have future event horizons characteristic of the de Sitter model. This means 
some events have causal futures that do not share any common events. In these dark energy scalar field models, some 
correlations are therefore unmeasurable, which destroys the observational meaning of the S-matrix. This indicates 
that it is not straightforward to bring superstring/M theory into consistency with dark energy models in which the 
expansion of the universe is accelerating 120 . 

At the time of writing, while there has been much work, it appears that the dark energy scalar field scenario still 
lacks a firm, high energy physics based, theoretical foundation. While this is a significant drawback, the recent flurry 
of activity prompted by developments in superstring/M and brane theories appears to hold significant promise for 
shedding light on dark energy. Whether this happens before the observations rule out or "confirm" dark energy is an 
intriguing question. 
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